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We study a gauge-invariant variational framework for the Yang-Mills vacuum wave functional. 
Our approach is built on gauge-averaged Gaussian trial functionals which substantially extend pre- 
viously used trial bases in the infrared by implementing a general low-momentum expansion for the 
vacuum field dispersion (which is taken to be analytic at zero momentum). When completed by 
the perturbative Yang-Mills dispersion at high momenta, this results in a significantly enlarged trial 
functional space which incorporates both dynamical mass generation and asymptotic freedom. After 
casting the dynamics associated with these wave functionals into an effective action for collections 
of soft vacuum-field orbits, the leading infrared improvements manifest themselves as four-gradient 
interactions. Those turn out to significantly lower the minimal vacuum energy density, thus in- 
dicating a clear overall improvement of the vacuum description. The dimensional transmutation 
' mechanism and the dynamically generated mass scale remain almost quantitatively robust, how- 

ever, which ensures that our prediction for the gluon condensate is consistent with standard values. 
Further results include a finite group velocity for the soft gluonic modes due to the higher-gradient 
' corrections and indications for a negative differential color resistance of the Yang-Mills vacuum. 

o 

; I. INTRODUCTION 

The variational approach was among the first theoretical methods developed to study nonperturbative ground-state 
properties of quantum systems [lj]. Soon after the inception of QCD, it was brought to bear on the physics of the gluon 
Oh| vacuum as well 0d^j j making increasingly use of its capabilities to treat nonperturbative problems analytically, at real 
time and at strong coupling (in contrast e.g. to semiclassical approximations). Indeed, variational techni ques remain 
applicable even in non-static and curved spacetimes @, H| as well as in nonequilibrium situations |l0||. Since 
variational problems are naturally formulated in the Schrodinger picture, furthermore, they often invoke quantum- 
mechanical intuition to shed new light on field-theoretic states, as encoded e.g. in the node number of their wave 
functionals [1, HH or in the action of transformation groups Hi. O n the other hand, field-theoretic applications have 
to implement the UV field modes and their dynamics exactly [ll| , they must deal with the standard renormalization 
issues in the Schrodinger representation [l2l [l3j , and they require a manageable way of calculating matrix elements by 
functionally integrating over physical intermediate states. For gauge theories, maintaining gauge invariance becomes 
an additional mandatory requirement which, in the Hamiltonian formulation, amounts to preserving Gauss' law. 
Traditionally, the latter is implemented in combination with (almost) complete gauge fixing to Coulomb gauge @, 
UH, although generally without full account of the gauge-group topology and the resulting Gribov copies. In addition, 
this gauge-fixed approach has to deal with a nonlocal Hamiltonian [161 ] and must take special precautions to avoid 
gauge-symmetry breaking approximations [l?} • 

Alternatively, the variational problem can be formulated in a (residual) gauge-invariant manner which, at least under 
certain approximations, remains analytically tractable (l7j . The mechanism of dimensional transmutation becomes 
manifestly gauge-invariant and particularly transparent from this perspective [l8[, and the dynamical mass gap (l9j 
emerges already from a minimal parametrization of the Gaussian "core" wave functional [l7| . The projection onto the 
gauge-singlet component of this core functional may furthermore generate an area law for spacelike Wilson loops and 
thereby account for linear quark confinement [l7l [Til . [20| , although experience from compact electrodynamics suggests 
that the vacuum wave functional will need non-minimal adjustments to describe the charged sector 21]. Another 
attractive feature of this approach is that it reexpresses the infrared dynamics in terms of gauge- invariant matrix fields 
(which subsume whole families of gauge- field orbits [22|) and that it preserves traceable links of heuristic value between 
these soft collective fields and the underlying Yang-Mills fields. Moreover, the gauge-invariant reformulation preserves 
the global, i.e. topological properties of the gauge fields. The contributions from all vacuum homotopy classes as 
well as instanton-mediated transitions between them, for example, are explicitly and transparently taken into account 
without recourse to the semiclassical approximation [51 Il7l. l22l l23j] . Besides such sets of (multi-) instanton and meron 
orbits, the dynamics contained in the vacuum wave functional sustains a rich variety of additional gauge-invariant 
infrared modes, including Faddeev-Niemi knots and other topological as well as nontopological solitons [22. l24j. 

Our main objective in the present article will be to improve the vacuum description provided by the gauge-invariant 
variational approach. To this end, we develop a trial functional family which accommodates a rather general mo- 
mentum distribution for the infrared vacuum modes while remaining analytically tractable. Due to the explicit 
representation of the ground state, structural aspects of the Yang-Mills dynamics and their impact on the vacuum 
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can then be understood in terms of the resulting gauge-field population. Moreover, our generalized trial basis will aid 
in the full exploration of the gauge-projected Gaussian functional space, e.g. by evaluating and analyzing relevant 
amplitudes, and thereby in the assessment of its limitations and further improvement potential. In the longer run, 
this type of analysis should reveal the extent to which Yang-Mills vacuum physics can be captured by the most general 
Gaussian core functionals, their subsequent gauge projection and the mean-field treatment of the resulting soft-mode 
dynamics. One may hope that these insights will eventually help to find a systematic and manageable approximation 
to the Yang-Mills vacuum wave functional beyond the Gaussian ansatz. 

After having set up our extended trial functional space, the next major task will be to establish the framework 
for calculating the associated vacuum amplitudes. This will enable us, in particular, to evaluate the vacuum energy 
density and its reduction in the generalized trial space, thus providing a quantitative measure for the improvement 
of the vacuum description. The analysis of the infrared mode distribution encoded in the energetically favored wave 
functional will further shed new light on the resulting vacuum physics. It will, in particular, reveal the physical 
significance of the variational parameters and clarify the impact of their optimized values on the soft-mode dynamics. 
We will furthermore examine the influence of the enlarged trial space on the phase structure and the order-disorder 
phase transition of the underlying infrared dynamics (which has, as the Yang- Mills dynamics, its origin in the imposi- 
tion of gauge invariance). In fact, the robustness of this phase transition ensures a reliable variational determination 
of the ground-state energy. Our analysis of the resulting vacuum properties will be complemented by evaluating 
the lowest-dimensional gluon condensate, which also provides a quantitative test of the dynamically generated mass 
scale and its relation to the underlying trace anomaly. Beyond the evaluation and analysis of specific amplitudes, 
finally, the extended variational framework will furnish a transparent theoretical laboratory well-suited to build up 
nonperturbative real-time intuition about the Yang- Mills vacuum [6l| . 

The paper is organized as follows: in Sec. |TT]we introduce our infrared-generalized, gauge-invariant trial functional 
family for the Yang-Mills ground state. In Sec. IHII we rewrite the associated vacuum overlap amplitude in terms of an 
effective bare action, summarize our strategy for calculating the variational bound on the Yang-Mills vacuum energy, 
and review the perturbative integration over the hard modes which results in a renormalized soft-mode action. In 
Sec. IIVI we develop the theoretical framework for calculating relevant amplitudes by integrating over the soft modes, 
and in Sec. [V] we derive explicit expressions for the 2- and 4-point functions. Their further evaluation is prepared 
in Sec. IVII where we set up the saddle-point expansion for the integral over the soft modes. We further derive the 
corresponding gap equation and find its solutions, which provide quantitative information on the phase structure of 
the soft-mode dynamics. In Sec. I VIII we calculate the various contributions to the vacuum energy density (both in the 
ordered and disordered phase) which we minimize in Sec. IVIIII We further compute the gluon condensate, discuss 
the physical implications of our results and suggest a few directions for future work. In Sec. IIX1 finally, we collect 
our principal results and our main conclusions. 

II. VACUUM WAVE FUNCTIONAL 

The success of variational estimates depends mostly on the choice of the underlying trial function(al) space. Ide- 
ally, it should be both comprehensive enough to accommodate the relevant physics and concise enough to remain 
analytically tractable and reasonably transparent. In the following sections we establish the basis of our subsequent 
analysis by introducing a gauge-invariant trial-functional family for the Yang-Mills ground state which is designed 
to compromise between these conflicting goals. Although our trial space expands the minimal one used in Ref. [I?} 
substantially, the ensuing variational analysis will turn out to require only a moderate increase in computational 
effort. 

Our basis implements Feynman's three requirements for trial wave functionals of gauge theories, i.e. the at least 
approximate calculability of matrix elements, the correct UV asymptotics and gauge invariance [TTj . The major 
methodological advance of Ref. 17] was to show that a trial space of manifestly gauge-invariant wave functionals, 
namely gauge-projected Gaussians, still permits an (approximately) analytically manageable variational analysis. In 
the following sections we first discuss the essential features of the minimal ansatz of Ref. [171 , and then generalize its 
momentum distribution for the infrared vacuum modes. 



A. Gauge invariance by projection 

Variational analyses are naturally set up in the Hamiltonian formulation of Yang-Mills theory and in the "co- 
ordinate" Schrodinger picture. This restricts gauge transformations to a fixed reference time, thereby effectively 
decoupling them from the dynamical time evolution. In the temporal (or Weyl) gauge Aq — which we adopt in 
the following, the residual gauge transformations are static, furthermore, and problems with negative-norm states 
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(as e.g. in covariant gauges) are avoided from the outset [62j. The gauge invariance of the vacuum wave functional 
is imposed by Gauss' law which acts as a subsidiary condition in Fock space. This is particularly crucial for varia- 
tional treatments because the Yang-Mills Hamiltonian becomes unique only when acting on physical states. Hence, 
as pointed out in Ref. [l7j . trying to minimize the vacuum expectation value of the Hamiltonian in the whole space 
of normalizable functionals could favor contributions from unphysical parts of the Hilbert space which may be almost 
arbitrarily enhanced by adding functionals of gauge-group generators to a given Hamiltonian. 

Starting from an approximate and hence typically gauge-dependent "core" functional ipo to be specified in Sec. Ill Bl 
(i/)q depends on the static gauge fields A (x), i.e. on half of the canonical variables), we impose gauge invariance by 
projecting on its gauge-singlet component, which amounts to averaging over the gauge group (26j |. The result is a 
trial vacuum wave functional of the form 
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where dfj, is the invariant Haar measure of the SU(7V C ) gauge group [63j . Q is the homotopy degree or winding number 
of the group element and 9 is the vacuum angle. The expression ([1} acquires the obligatory 9 phase under large 

(i.e. topologically nontrivial) gauge transformations and renders Gauss' law manifest. 

Instead of dealing with the gauge-invariant trial functional (JlJ, one may alternatively evaluate the energy density 
of gauge- dependent Gaussian wave functionals and then correct for the lack of gauge invariance before minimiziation. 
This approach was pursued in Ref. [27] where the spurious kinetic energy due to gauge rotations was subtracted 
by Thouless-Valatin projection (as originally developed for the treatment of deformed nuclei). The main advantages 
of the approach based on wave functionals of the type ([!} are that gauge invariance of subsequent calculations and 
approximations is maintained exactly, and that it applies universally to all matrix elements. 



B. Gaussian core wave functional 



In order to determine the gauge-invariant trial functionals ([T} and the related amplitudes completely, it remains 
to adopt a core functional family ipQ. As in any variational calculation, this choice generally remains an uncontrolled 
approximation without systematic improvement strategy. In order to motivate our choice for ip and to discuss the 
underlying physics, we first recall that the vacuum wave functional cannot have nodes, i.e. that ipo [A] > for all A j4j]. 
(This provides an example for how quantum-mechanical insights remain applicable to field theory in the Schrodinger 
representation as long as its infinitely many degrees of freedom do not generate qualitatively new effects.) Hence one 
may write without loss of generality 
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where Af is a generally infinite normalization constant and $ is a real and typically nonlocal functional of the gauge 
field. In order to find a viable approximation for $, we represent it as a functional power series in A and determine 
physically reasonable and analytically tractable truncations. The first two terms of this series are generally discarded: 
constant terms can be absorbed into the normalization constant Af, while the term linear in the gauge field (by itself) 
corresponds to a coherent vacuum state which is known to be unstable [28| . 

The next term is quadratic in A and plays several crucial roles. The first originates from an ambiguity in $ which is 
due to the invariance of the Haar measure in Eq. (JTJ) under right (or left) multiplication by any given group element. 
As a consequence, infinitely many choices for $ [A] lead up to an unphysical redefinition of the normalization constant 
to the same ^>o [201 ] - This ambiguity can be removed, however, by prescribing the quadratic term of $. Furthermore, 
this term leads to a product of Gaussian vacuum wave functionals of the Abelian U(l) gauge theory [64| which 
asymptotic freedom renders exact in the ultraviolet. Finally, and from the practical perspective most importantly, 
the Gaussian functional ipo resulting from a quadratic term can be integrated over A analytically, while higher-order 
contributions may at best be treated perturbatively [3fJ. Hence one generally truncates the series for $ after the 
quadratic term, which leads to the "squeezed" approximation for the core functional, i.e. to the Gaussian 
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with the normalization factor A/q 1 = [det (G/2)] 1 ^ 4 and a real kernel or "covariance" G -1 which satisfies a normal- 
izability condition (cf. Eq. ©). Eq. © appears to be the "richest" core functional family whose matrix elements 
can be dealt with by the currently available analytical methods of field theory. (Adding c number sources to the 
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variable A would still allow for an analytical treatment and generate finite vacuum expectation values for A in gauge- 
fixed approaches, while such local "condensates" (except of the time component) are not sustained in our residually 
gauge- invariant vacuum and would be rendered mute by the gauge projection in Eq. (JTJ.) 

As already mentioned, the core functional © represents an infinite product of Gaussians, one for each Fourier mode 
of the gauge field with momentum k. Hence the components of the Fourier-transformed covariance 

< $ : = G t' ah = 9? (*) % + 9? (k) kk (4) 

(where fcj := ki/k with k := \fkjki) turn after diagonalization into mode frequencies or energies. The squeezed state 
thus generalizes the ground state of the quantum-mechanical harmonic oscillator. In our context, it corresponds to a 
vacuum consisting of color-singlet gauge-field pairs which may be regarded as the protostate of a glueball condensate 
(Sjj . (This should be contrasted to the "condensation" of single gluons in a coherent state.) 

Since Gaussian functionals transform nontrivially under non-Abelian gauge groups, the gauge averaging in Eq. ([1]) 
is an integral part of our approximation to the physical vacuum state. Gauge-fixed variational analyses in Coulomb 
gauge (for references see e.g. 0, H, 1(1 JliL Jl5l | ) are generally based on Gaussian functionals as well, however. They 
were found to generate a mass gap [2|, Il4j and, when multiplied by the inverse square root of the Faddeev-Popov 
determinant, an approximately linearly rising confinement potential [l5j . 



C. General properties and UV asymptotics of the covariance 

We are now going to specify the properties of the covariance ((4]) which characterizes the members of our core trial 
functional family and which contains the variational parameters whose values will be adapted below to optimally 
approximate the Yang-Mills vacuum wave functional. Translational invariance implies G~ 1 ' ab (x,y) — G^- 1 '" (x — y) 
and was already anticipated in Eq. (|3]). Without loss of generality at the perturbative level (and beyond, since the 
integration over the gauge group in Eq. 0} averages out longitudinal contributions), we will further specialize to a 
purely transverse covariance with g^ = [17| , i.e. 

G^ ab (k)=S tJ G- 1 ' ab (k), (5) 

which allows for a direct comparison with the results of Ref. [l7| and should be sufficient for our explorative purposes. 
(The impact of longitudinal contributions cx g% b was discussed in Refs. [l8|, HII, and a specific prescription for g% b was 
shown to reproduce the 1-loop Yang-Mills /3-function. Due to the ambiguity of the core functional mentioned in Sec. 
IIIB1 a longitudinal term in the covariance can always be removed from $ without changing the physical part of the 
vacuum wave functional, although generally in exchange for terms containing higher powers of the gauge field (20|.) 

In order to motivate the color structure of our covariance, we note that only the homogeneous part of the gauge 
transformations keeps the exponent $ of the core functional bilinear in A. Gauge transformations which vary little 
over distances for which G _1 (a? — y) has appreciable support therefore leave the Gaussian part of the exponent with 
a covariance of the form [TtJ 

G~ l ' ah (k) = d^G- 1 (k) (6) 

approximately invariant. This will hold, in particular, for those gauge group elements which determine our soft-mode 
dynamics (see below). Since the same color structure is appropriate for the hard gauge- field modes, which we will 
treat perturbatively in the small bare coupling c/b, we adopt Eq. (|6]) for the remainder of this paper. As a consequence, 
the core functionals ipo become invariant both under global U(N C ) transformations and under N% — 1 copies of the 
U(l) gauge group. 

Combining the spacial ([S]) and color ([5]) structures, our covariance assumes the form 

G^ ab (f, y) = S^G- 1 (x-y) = 5 ab 5 l3 J ^^^G^ (k) (7) 

which is symmetric in each of the 3 "index" pairs. The covariance is further restricted by the requirement that wave 
functionals of physical states have to be normalizable. Since the norm of the functional ([3]) involves an integral over 
the Fourier modes A (fc), this implies that the integrand ipQ (k) ipo (k) has to damp large gauge fields A (k) for all k. 
The corresponding localization in field space is implemented by demanding 

G- 1 (k) > (8) 
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which ensures vacuum stability and a positive energy spectrum of the associated quantum field theory. The condition 
([5]) will limit the domain of the variational parameters to be introduced below. 

As noted by Feynman, a further mandatory requirement is that all trial functionals reproduce the asymptotically 
free gauge dynamics for k — > oo exactly . This prevents the infinitely many ultraviolet modes (which are irrelevant 
for the vacuum physics) from artificially dominating the soft-mode energy density through their IR-modc couplings. 
In order to implement the correct UV behavior, we factorize the unprojected core functional ([3]) as 



4 G) A 



4 G<) - 1 



(9) 



by splitting the k integration domain in the exponentials into soft/hard regions with momenta k ^ 
separation scale /i, i.e. 



55 /i relative to a 



(%) 



exp 



d 3 k 



Tk 2 )A^ (fc)G^ 1 (k)A^. (k) 



(10) 



This allows to incorporate the asymptotic freedom of Yang-Mills theory (i.e. the Gaussian UV fixed point) by requiring 
that G approaches the non-interacting, massless static vector field propagator Go (k) = 1/k when k — > oo. As long as 
[i 3> Aym where Aym is the Yang-Mills scale, a natural approximation is therefore [65| 



GZ, 1 (k) = k 



(11) 



which we adopt in the following [66|. More specifically, the value of /i must be chosen large enough for perturbative 
corrections from the hard modes, or equivalently the renormalization-group (RG) improved running coupling a (p) = 
g 2 (/i) / (47r), to remain small. (RG-improved perturbative corrections to the leading k — > oo behavior will enter 
when integrating out the hard modes perturbatively, cf. Sec. IVII Al and App. IA 21 ) Since /i will be treated as a 
variational parameter, this has to be checked a posteriori, i.e. for the value /x* which turns out to minimize the 
vacuum energy. For the infrared momenta fc < /i, finally, the nonperturbative Yang-Mills dynamics is expected to 
induce a qualitatively different covariance G^ 1 (fc) which will be the subject of the following section. 



D. Generalized IR (soft-mode) covariance 

Due to the logarithmically slow running of the Yang-Mills coupling in the ultraviolet, the high-momentum behavior 
of our trial functional family is rather accurately reproduced by the hard-mode covariance (fTTj) . Hence only the IR 
covariance G<* (fc), which encodes the more complex and less understood nonperturbative vacuum physics, remains to 
be determined variationally. To this end, we have to implement a parametrization for G^ 1 which is sufficiently "rich" 
to accommodate the relevant physics without impeding an analytically tractable variational analysis. The minimal 
choice 
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was adopted in the pioneering study [171 ] and shown to generate a mass gap when the separation scale /j, is simulta- 
neously utilized as a variational parameter. While it therefore provides an efficient starting point for describing the 
Yang-Mills vacuum, one may also worry about too much bias since it leaves only one characteristic vacuum mass scale 
to be determined by energy minimization. This involves the risk of essentially "building in" the mass gap without 
gaining much further insight into the underlying vacuum structure. 

In the following, we will therefore rely on a more comprehensive parametrization of G^ 1 which better accommodates 
the k < \i mode dynamics while still allowing for a variational analysis without the need for solving a functional 
differential equation. It is based on the under reasonable analyticity assumptions general gradient expansion [22[ 



GZ l {x-y) 



d 2 

1 i x 



'd 2 x 

5 { y 



6% (x-y) 



(13) 



which can be efficiently truncated to yield a manageable trial basis for the k 2 <§; /i 2 soft-mode physics. Besides /i, the 
variational parameter space then contains the IR gluon mass m g > and a few of the low-momentum constants c n 
which characterize dispersive properties of the vacuum (cf. Sec. IVIII CI) . (The parameters c„ should be considered as 
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renormalized at fi since they do not receive UV contributions from integrating over the k > fi modes.) The derivatives 
in Eq. (fT5|) act on the regularized delta function 
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which encodes the slow variation ||9A< || / || A< | < /i of the soft modes and ensures that the higher-order terms of the 
above expansion are parametrically suppressed. It further restricts the support of G^ 1 : 
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G< x (£) ~ for |f| > - 



(15) 



Together with the hard-mode covariance (jllj). Eq. (|13j) provides a rather general paramctrization of the core functional 
(|3"|) which includes the minimal one-parameter ansatz (1121) for c n = and m g = //. 

(Note, incidentally, that the perturbative covariance (fTTj) has a branch point at fc 2 = and thus cannot be directly 
expanded as in Eq. (H"3l) . Note further that the expression (H"3l) can be extended beyond the limited spacial tensor 
structure ([5]) by writing 



G<^- (x - y) = G^l; 1 (x - y) + G^ X (f - y) 



(16) 



with (X = L, T) 
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kikj/k 2 . This decouples the longitudinal and transverse contributions 



0. 



(19) 



For m 
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the more general expression (fTS)) simplifies to Eq. (fT3"|). In the following analysis 



we adopt the latter, in order not to compromise transparency by additional parameters and to allow for a direct 
comparison with the results of Ref. [l7|.1 
In momentum space Eq. (|13[) becomes 
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which renders the relation between the parameters c„ and the dispersion properties of the IR quasigluons (cf. Sec. 
IVIII CP more explicit. It further implies that the effective IR gluon propagator 
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resulting from G< (fc) G< (fc) = 1, is analytic at fc 2 = 0. (The finiteness of G< at zero momentum, G< (0) = m~ 1 , 
is reminiscent of lattice results for the gluon propagator in Landau gauge [33j.) The adjustable parameters m g , fi 
and Cn are to be determined variationally (cf. Sec. IVIII[) . Their physically sensible domain is restricted by several 
constraints, however. Indeed, as a consequence of the normalizability condition ([SJ the low-momentum constants 
must satisfy the bounds 



C\ < 1, C 2 > -1, 



(22) 
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(for TO g > 0). The value of the IR gluon mass m g > is further constrained by requiring continuity of G _1 (k) 
at the matching point k — /j, between soft and hard momenta. This fixes m g as a function of the other variational 
parameters. When approximating the series (|20| by the truncation c n >2 = 0, for example, one has 

TO g 0»,ci) = — . (23) 
i — a 

Note that the requirement of a non- negative IR gluon mass then restricts the c\ domain to C\ < 1, in agreement 
with the bound (|2"2"|) from normalizability. The singularity of m g ([i, c\) for c\ — > 1 reflects the onset of the vacuum 
instability and is inherited by G^ 1 (k), cf. Eq. (|20[) . 

In addition to the lowest-order approximation (|12p , our expansion (| 13[) encompasses another previously used ansatz 
for the IR covariance, namely the inverse of the non-interacting massive vector propagator, 

G 1 (k; f i) = Vfc 2 +M 2 k = A' 

which obviously corresponds to m g = /i, = — 1/2, c^' -1 = —1/8 etc.. Its truncation to c n >2 = was adopted as 
the basis of a gauge-invariant saddle-point expansion for the IR amplitudes and found to reproduce the soft-mode 
action (cf. Sec. IIIII) at its saddle points within a few percent accuracy [22j. (Lattice simulations similarly found the 
Landau-gauge gluon propagator to behave like a massive vector propagator at intermediate momenta |34j|.) 

In Fig. [Tjwe compare the expressions (fT2"j) and (pM|) for the IR covariance to Eq. (|20p with ci — ±0.15 (the positive 
value will turn out to minimize the vacuum energy, cf. Sec. IVIII A[) and c n >2 = 0. A more detailed discussion of the 
physics encoded in the covariance (|20p will be postponed until the variationally optimized values of the low-momentum 
constants are found (cf. Sec. IVIII C|) . 

In summary, the Gaussian ansatz ([3]) for the core functional encodes the correct ultraviolet asymptotics, can be 
partly motivated in the infrared and is amenable to existing analytical methods of field theory. As such, it provides 
the best currently available, analytically manageable and gauge-invariant approximation to the Yang-Mills vacuum 
wave functional. 
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III. ENERGY AND SOFT-MODE DYNAMICS OF THE TRIAL FUNCTIONAL FAMILY 



A. Calculational strategy 



Variational analyses a la Rayleigh-Ritz amount to minimizing the expectation value of a given Hamiltonian in a 
suitable space of trial states. In our context, the energy density of the trial vacuum is 
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where % is the Yang-Mills Hamiltonian (cf. Eq. (113110 in temporal gauge and the generalized vacuum wave 
functional family introduced in Sec. [TTJ This expectation value is most efficiently evaluated by means of a generating 
functional, as summarized in App. [Aj In the present section we prepare for the calculation of Eq. (|25|) by sketching 
our basic strategy (which straightforwardly generalizes to more general matrix elements, cf. App. lAl . 

After inserting the gauge-projected vacuum wave functional ([T]) into Eq. (|25|) and interchanging the order of the 
integration over fields and group elements, the gauge invariance of the A integral allows to factor out one gauge group 
volume. Eq. (|2"5)) can be then be rewritten as 
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(where DU is the functional SU(A r c ) measure as defined in Eq. ([TJ). After evaluating the functional derivatives 
contained in H, the Gaussian integration over A can be performed exactly (cf. App. [A] for more details), resulting in 
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(27) 
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where wc introduced the abbreviation 

/ l...A...E...E)))exp{-T h [[/]} 
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The above expression defines an effective bare action r b [U] which describes dynamical correlations originating from 
the gauge projection of the functional tpo. More specifically, from the normalization (((!))) = 1 it follows that 



T b [U] = — In / DA ipQ 
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This gauge-invariant, 3-dimensional Euclidean action would become U independent if ipo were gauge invariant by 
itself. Hence r b [U] gathers all those gauge- field contributions to the generating functional whose approximate vacua 
ipo at t — ±oo differ by a relative gauge orientation U. The variable U thus represents the contributions from a 
specifically weighted ensemble of all gluon field orbits to the vacuum overlap and is gauge invariant by construction. 

After splitting U (x) = U < (x) E7> (x) with £/> (x) = exp (—ig(j) a (x) X a /2) into hard- and soft-mode contributions 
and integrating over the hard UV modes cj> perturbatively (as done in Ref. |17| which used the same UV covariance 
(|TT|) . cf. App. 12}, one arrives at 



(H) 



J BU < J (((H))) exp {-r b [0, t/<]} 
r£>E/< f£^exp{-r b [^E/<]} ' 



After further defining 



((O)) exp {-r< [[/<]} 



D</>(((0)))exp{-T b [</>,£/<]} 



(30) 



(31) 



(32) 



we end up with a reformulation of the matrix element (1251) in terms of the dynamics for the t/< field, i.e. for the 
low-momentum components of the integration variable originating from the gauge projection of the vacuum functional 



and in particular the effective soft-mode action 

r<[C7<] --In/ £>0exp{-r b [</>,£/<]}• 



(H) 



\DU < ((U)) exp {-r< [[/<]} 
/ DU < exp {-r< [[/<]} ' 



(33) 



Hence the calculation of (H) boils down to an integral over U < in which the "reduced" (i.e., fixed [/<) matrix element 
((H)) is weighted by a Boltzmann factor containing just the soft-mode action T< [£/<]. This action is governed by 
our IR expansion (fl3|) of the covariance G^ 1 which produces interactions to be determined in Sec. IIIIBI Since ((H)) 
consists of nonlocal functionals of the soft modes U < , furthermore, the evaluation of Eq. (|3"3"]l amounts to calculating 
(equal-time) soft-mode correlation functions. The necessary framework for such calculations will be set up in Sec. IIVI 



B. Soft-mode action 



In this section we derive explicit expressions for the low-momentum dynamics (|32[) induced by the wave functional 
(P) with the Gaussian core © and the generalized covariance (fTTl) . (|2"01) . After specializing the expression (|2"9")l for 
the bare action to the core functional ©, the integral over the gauge fields can be carried out exactly (cf. App. [A"l) . 
The result is the action of a 3-dimensional, bilocal nonlinear sigma model: 

Tb [U] = ^ J d 3 x J d 3 yL<f (x) V ab (x - y) h\ (y) . (34) 

(Above we have omitted a term of higher order in the small bare coupling g^, cf. App. [A]) The U fields enter r b 
both in terms of the left-invariant SU(A C ) Maurer-Cartan forms 



Li(2) = CTt (5)fttT(jg)=:£f (x) ^ 



(35) 
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(with real components Lf and the SU(7V C ) Gell-Mann matrices A a ) and through higher-order corrections to the bilocal 
operator 



ab 



(G + G l 



ab 



_ G -l S ab 



(36) 



where G u = G ab (x - y) W [x] t a U (x) ® U{y)t b W (y) and t a = A Q /2. The above reformulation of the vacuum 
functional overlap can alternatively be obtained by a saddle-point evaluation [18( of the Gaussian integral in Eq. (|29l) . 

According to the strategy outlined in Sec. IIII A[ the calculation of infrared amplitudes involves an effective action 
(f32|) which only retains soft field modes explicitly. After factorizing U = U < U > , the hard modes U_> with k 2 > fj, 2 can 
be integrated out of Eq. (f3"2"| perturbatively (cf. App. [AJ as long as the renormalized coupling [35| 



f ^ i ffb^c , A uv 
9 W = 9b + - — — In ■ 



(27T) 



/' 



(37) 



(for G> (fc) = fc _1 ) stays sufficiently small. Note that Eq. (|37l) exhibits asymptotic freedom and differs from the 
one-loop Yang-Mills coupling only by a small correction factor 1/11 which could e.g. be generated by an anisotropic 
component for Gz} [l8l . |32| . The one- loop integration over the high- momentum modes was found to be reliable down 
to /i ~ 1.3 GeV [35j . In this range, the resulting renormalized soft-mode action 



r< [u<] = 



i 



V (m) 



d 3 x / d 3 yL a < Jx)G- 1 (x-ff)L a < Aff) 



(38) 



(cf. Eq. (|A48[1 ) is obtained from Eq. (1341) by replacing the bare coupling gi, with the running coupling g (//). One may 
further simplify the action ([3"5| by exploiting the strongly local support of the IR covariance (fl"3"|) (compared to the 
minimal soft-mode wavelength fj,^ 1 ) which manifests itself in the regularized delta function (1141) and its derivatives. 
Indeed, this allows to apply the unitarity constraint for the U < approximately at neighboring points, i.e. 



G- 1 (x - y) [4 (x) C/< (y) ~ G' 1 (x - y) , 



(39) 



which could be systematically improved by Taylor expansion of the slowly varying soft modes. To leading order, one 
may therefore replace Eq. (f3"5)) by the 2U contribution 



r 



(2U) 



[U<} = 



1 



2ff 2 (M) 



d 3 ytr [d i U < (x) G' 1 (x - y) (y)} 



(40) 



Inserting our full expansion (|13[) of the generalized IR covariance into the low-momentum dynamics (|38p. on the other 
hand, yields the complete soft-mode action in the form 



T< [[/<] = J d 3 z [£ <i0 (z) + £ <>C1 (z) + £ <iC2 (z) + ...] 



(41) 



where the £< jC „ (z) are (quasi-) local Lagrangians. (A saddle-point expansion of this dynamics, with the c„ determined 
by the massive vector covariance (T2"4"]) . was employed in Ref. 22] to identify gluonic IR degrees of freedom.) The 
leading-order, two-derivative Lagrangian 



£< '° (f) = M (f) 9lU< ^ U < {g) dlU< = ^ { d i U < d * U <} 



(42) 



(where we used U' < U < = 1 and neglected a surface term) is just the standard nonlinear a model. The c„ ^ 
corrections (with the low-momentum constants c„ restricted by the requirements of a positive static action and a 
bounded vacuum energy, cf. Sec. Ill Dp generate 2 (n + l)-derivative interactions which are, relative to the leading 
term (|4"2"]) . suppressed by n powers of k 2 /[i 2 : 



(*) = ^ktr {c/t {g) diU< [g] d m ^.^t (f) ^ I 



(43) 



The action (jUJ), based on the Lagrangians (|4^)) . (|43l) . exhibits the dynamics generated both by the generalized IR 
covariance (|13[) and by gauge projection. It encodes, in particular, information on the vacuum topology which shows 
up in the form of instanton (3(| (through the Atiyah-Manton holonomy 37 1), meron (38|, monopole [39[, Fadeev- 
Niemi-Cho knot [4(J |4l| etc. contributions and can be made explicit 
reduce the c n ^ interaction terms to their bilinear (i.e. 2U) parts 



C 



(21/) 



2ffV 



The approximation ([391 can be used to 
tr[diUl (zOd 2 ™^ (£)]}, (44) 
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furthermore, which will turn out to generate the dominant contributions to the vacuum energy density and other 
matrix elements. 

We close this section by recalling that the soft-mode dynamics (|41j) follows uniquely from the adopted vacuum wave 
functional and preserves traceable links between the ?7< fields and the underlying gauge fields |22|. Nevertheless, 
it remains reasonably transparent and allows for efficient and controlled truncations resulting e.g. in the analytical 
expressions for the vacuum energy to be derived below. These benefits originate in large part from reexpressing the 
dynamics in terms of the gauge-invariant fields U which gather collective contributions from whole gauge-field orbits 
instead of dealing with each gauge field individually [221 ] . 



IV. SOFT-MODE CORRELATION FUNCTIONS 



As outlined in Sec. IIII Al the calculation of the vacuum energy density requires the evaluation of the reduced matrix 
element ((H)) which contains correlations functions of the soft- mode fields [/<. In the following section we set up our 
framework for evaluating these correlators. 



A. Generating functional 



Soft- mode correlation functions, as they appear in the integrand ((H)) of the vacuum expectation value (|33[) and 
in the gap equation to be derived below, are efficiently calculated by means of the generating functional 



Z [j,J + ] = J DU < cxp -r< [[/<] - J d 3 ztr [jU{ + j+E/< 



(45) 



where are matrix sources. In order to prepare for the (approximate) integration over the soft modes in Eq. (|45[) . 
we first release the unitarity constraint on the U < fields in the usual manner by inserting a delta functional, i.e. 



Z = J DVS [yV-l] exp -r< [V] - J d 3 ztr{jVl + j f V} 

= J DE J DVexp -r< [V] - T s [V,S] - J d i ztr{jV' i + jV} 



(46) 
(47) 



where the Hermitian matrix fields E (x) act as Lagrange multipliers (with a normalization chosen for later convenience) 
and where Fs [V, E] contains the interactions between the E and V fields, 



T E [V, E] = J d 3 x£z (z) = g§ J d 3 xtr [E (V*V - l)] 



(48) 



The integral over the unconstrained complex matrices V has a linear Euclidean measure DV. Hence, after pulling out 
the subleading (cf. Sec. IIIIBI) . non-Gaussian AU interactions as functional derivatives with respect to the sources, it 
can be performed analytically. This results in 



Z[j,J ] ] = J DEW exp j- j d 3 z[C < . + C s +C 

( r {W) 
\ 1 ~<,ci 



DYjDV exp 

y(«0 



5_ _S_ 



d 3 z 



(4C7) 
<,c 2 



£ <)C2 + ... + tr {jV*+jiV}] 

d 3 z(c^ +tr{jV^ +jV}) 



exp 



where we have gathered the part of the generating functional which originates from the bilinear Lagrangian 
C,W) : = C <fl + C { Z\ + £ { < U d + + £s =: tr {V^AV} ^tr {E} 



with the kernel 



2.9 2 



A(x-y; {<;,} ,E) := ^ + + + ... - E ) 6% - //) 



(49) 
(50) 
(51) 

(52) 

(53) 
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(which defines the inverse soft-mode propagator) as 

[j,j1]=[ DEW exp - J d?z(c^ +tr{jV^ +jV}) 
All contributions from the AU vertices, on the other hand, are collected in the functional potential 

V<«0 [u\u] =ex P 



'<,Cl ' J ~'<,C2 



The Gaussian integral over the V fields in Eq. (|54|) can then be performed analytically, with the result 

[j,jt] 



Z( 2C/ ) [0,0] 



£>Eexp 



d 3 x / d 3 2/tr {jt (f) A" 1 (f - y; {c„} , E) j (y)} 



(54) 



(55) 



(56) 



The whole E dependence of the correlators, as well as the c n dependence originating from the 2U interactions, is now 
concentrated in the gauge-invariant soft-mode propagator A -1 which we will analyze in Sec. IIV CI The perturbative 
treatment of the AU vertices for |c„| -C 1 starts from the expansion 



S_ _5_ 

Sj ' <5j t 



exp 
1 - 



d 3 AW) (S_J_ 



of the functional potential which contains the c n dependence induced by the AU interactions. 



(57) 
(58) 



B. Mean- field approximation and vacuum phases 



After having rewritten the generating functional as an integral over the auxiliary field E, we evaluate the latter in 
the saddle-point or mean- field approximation (cf. Appendix I A 31) . When the integrand in Eq. (|47[) is expressed as a 
Boltzmann factor 



exp 



{-f [£]} 



Wexp{-r<[y]-r s [v,s]} 



(59) 



(in the present section we are interested only in the vacuum overlap amplitude and therefore set all sources to zero), 
the saddle-point equation takes the form 



= g|expf[E] J DV [Ft (x)V(x)-l] exp{-r< [V] - r E [V^,E]} = 



(60) 
(61) 



which ensures that its solutions extremize the effective action T. To leading order in the saddle-point expansion, 
integrals over E are then approximated by their integrands where E is replaced by a solution E^™^ of Eq. (|6T|) . This 
turns Eq. (|ST|) into 



(Eft (x)U{x)) = l, 



(62) 



in particular, and thereby restores the unitarity constraint at the mean-field level. In a homogeneous vacuum one 
expects the solution of Eq. (|62|) to be a constant field. The interactions (|48|) then turn into a mass term for V which 
triggers dimensional transmutation (as in O (n) models) and generates a mass gap. Hence Eq. (|62]) plays the role of 
a gap equation. Since for E > the SUi(A r c )x S\Jr(N c ) symmetry of the soft-mode action (f4T|) is unbroken, the 
mean-field solution should further be proportional to the unit matrix, i.e. 



E( mf ) = E x 1, 



(63) 



where E is a real constant. 
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In order to exhibit the vacuum phase structure encoded in the wave functional (JlJ more fully, an analogous mean- 
field treatment of the V integral in Eq. ([4"Tf (which reproduces its Gaussian part exactly) turns out to be useful as 
well. Defining the associated effective action (again for vanishing sources) by 

exp{-f [V]} := J DEexp{-r<[Vl-r E [V;E]} (64) 

one obtains the corresponding saddle-point equation 

-expT[y] / DY, TTTTJ-a cxp{— r< [V] - Ts [V, E]} (65) 



SV^(x) ^ J 5Vi(x) 



-> ^| expr [V] [EV] exp {~r s [V, E] } = 0. (66) 

In the second line we have saturated the E integral with the constant (i.e. vacuum) mean fields V and E. This implies 
r< [V~\ = 0, in particular, since T< only contains derivative interactions. Hence Eq. (|66l) reduces to 

(EC/(f))=0. (67) 

(The SUl(A^ c )x S\Jr(N c ) symmetry of T (under which U,V transform as V — > RVL^) turns a nonzero mean field 
V into a continuous family of degenerate saddle points. The associated Goldstone zero-mode contributions are not 
suppressed by the Gaussian weight and have to be integrated exactly.) 

The solutions of the saddle-point equations (|62p and (|57]) characterize the vacuum phases of the effective a model 
(T4"Tj) . The gap equation (|^|) determines the auxiliary mean field E and will be solved in Sec. EH The saddle-point 
equation (|6T[) shows that the vacuum can exist in two phases, as expected on general grounds and confirmed by 



lattice simulations [42j, |43| and the e-expansion [44j. More specifically, with increasing "analog temperature" </ 2 (/i) or 
decreasing fi one expects the vacuum to pass through an order-disorder phase transition, as in the analogous statistical 
spin model. In the ordered low-temperature (i.e. weakly coupled) phase one has 

(U) ? 0, (E) = 0, (68) 

i.e. the SU l(N c )x S\Jr(N c ) symmetry is broken to its diagonal subgroup and — 1 massless Goldstone bosons are 
generated. In the disordered high-temperature (strong-coupling) phase with 

(U) = 0, (E) jt 0, (69) 



on the other hand, the symmetry of the action is restored [6Jj. The disorder-order transition occurs when (E) reaches 



zero in the disordered phase. As usual, it is the result of a competition between the ordering tendency of the energy 
and the disordering propensity of the entropy. Our above arguments imply, furthermore, that the qualitative phase 
structure is independent of the detailed interactions (1431) as long as the soft-mode dynamics contains only derivative 
interactions. (This is guaranteed by the unitarity of the U < fields.) As a consequence of T < [V] = 0, in particular, the 
mean-field equations and hence the qualitative phase structure have no explicit dependence on the low-momentum 
constants c„. 



C. Soft-mode propagator 

The static soft-mode propagator A -1 determines the Gaussian part (|56[) of the generating functional and thereby 
encapsulates most of the impact of the higher-derivative interactions on vacuum structure and amplitudes. From the 
Fourier transform 

a » " 2* / H 9 ' k3) - 7 k ' + 7 k °-- + £ ) e ' E '*~ i> (70) 

of its inverse (|53p. which we have specialized to the constant mean field (|63|) . the propagator is obtained as the solution 
of 

J d 3 zA- 1 (x-z)A(z-y) = S a < (x-y) (71) 
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one therefore has 



with 



where is the regularized delta function (fT4|) . In terms of the dimensionless variables 



*=-, e^— (72) 

A- 1 (x-y) = [ AjL e i^(S-y) A -t (k) (73) 
J (27r) 

x 2gV 0(l-* 2 ) 

1 J m ff k 2 + £ 2 + A4 2 (k 2 ) ' 1 j 

The Fourier representation of the soft-mode propagator can be simplified by performing the angular integrals analyt- 
ically, 

a-^x-v)--^ 1 t dn K8 ™M x -y\) (75) 

[X y) ~ ^m g \x-y\J d \^ + e+MH^Y {7b > 

which becomes useful for numerical implementations and demonstrates that the isotropic regularization (i.e. k < 1) 
preserves 

A- 1 (x-y) = A- 1 ($-x). (76) 
Equation (|74|) reveals that the higher-derivative interactions generate a momentum-dependent selfenergy 

M 2 (k 2 ) =-c iK 4 + c 2 k 6 -c 3 k 8 + ... (77) 

for the soft modes. Hence the propagator (|74|) may be regarded as the resummation of a static self-energy (i.e. M 2 ) 
insertion. This geometric series converges for \M 2 (k 2 ) / (k 2 + £ 2 )| < 1. Expecting the c\ contribution to dominate 
at small momenta k <C fi and therefore truncating to c„>2 = (see below), this implies 

k 2 -4- F 2 

\cl\ (78) 

K 

For k £ [0, 1] and ^ > 0, the above inequality is satisfied as long as |ci| < 1. This does not additionally constrain the 
parameter space, however, since c\ < 1 guarantees the normalizability of the vacuum wave functional (cf. Eq. (|22jl ) 
and |ci| <C 1 ensures the validity of our perturbative O (ci) treatment (1551) of the residual AU interactions. 
In order to analyze the singularity structure of the soft-mode propagator, we rewrite Eq. (|74[) for c„>2 = as 



[) m g K 2 {l- Cl K 2 ) + e ,„,, x i . UrC 2 1 ■ ;! ' 



which reveals two poles at the momenta 



«i,a (C,ci) = — (l±Vl + 4 Cl e 2 j. (80) 

The pole position k 2 (with the positive sign in Eq. (|80)) ) diverges for c\ — > but decreases monotonically for c — » 1 
to reach k 2 ci = 1) > 1. Hence the pole at k\ lies outside of the integration range k £ [0, 1] for all c\ £ [— oo, 1] 
(and even for £ — > 0). The second pole at k 2 . lies inside the integration range for c\ > —oo and £ = 0. For c\ = 0, in 
particular, it turns into the standard infrared pole of the c„ = propagator, showing that the momentum-dependent 
selfenergy does not create singularities beyond the £ — > pole of the uncorrected soft-mode propagator. 



V. EVALUATION OF THE 2- AND 4-POINT CORRELATORS FOR c n > 2 = 



Having established our calculational framework for the generating functional, we now derive explicit expressions 
for the soft modes' two- and four-point functions. The results will underlie our subsequent evaluation of the matrix 
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elements \ U< (x) U < (x)j in the gap equation (j52"j) and (£" i (x) L< ^ (y)) in the chromo-electric contribution (|142[) 

to the vacuum energy density. 

As stated previously, among the contributions from the higher-derivative interactions (|43l) those associated with c\ 
are expected to dominate at small momenta k 2 <C /j 2 . This expectation is further supported by evidence from Ref. 
(22j where relevant infrared physics was found to be captured by the truncation of the covariance expansion (|13p to 
c n >2 = 0. (More specifically, the first correction term £<. Cl turned out to reproduce the full (nonlocal) IR action (|3"5)) 
at its saddle points with better than 10% accuracy.) Furthermore, the c„>2 corrections (although straightforward 
to implement) obscure the (e.g. graphical) analysis of the vacuum energy density and the interpretation of other 
amplitudes. Hence we will restrict the remainder of our investigation to the c\ corrections, associated with the 
four-gradient interactions 



ClTOg 

,(2U) 



tr 



Ul{z)d i U < {z)d 2 



d l u{ (z) U. 



£^(z)+4 4 ^(z) 



(81) 
(82) 



(where derivatives are implied to act on the nearest field only) which we have split as in Sec. 1111 B I into the dominant 
2U interaction term £<^2, given by Eq. (|4"4"|) for n = 1, and the 4C7 contribution 



4 4 2 (*) 



c\m g 



tr 



U\d i U < 



2d l d k u\d k u < + d l u] i d 2 u. 



(83) 



As discussed in Sec. I IV A[ we will evaluate the 2XJ contributions to the relevant n-point functions exactly and treat 
the residual contributions from the ill vertices perturbatively [o8| to O (ci). The latter has to be justified a posteriori, 
by showing that the variational results favor small enough coupling values |c*| <C 1 (cf. Sec. IVHI A|) . 



A. The 2-point function 



In this section we evaluate the 2-point soft-mode correlator (for q>2 = 0) which simultaneously provides the basis 
for the calculation of the 4-point correlator in the next section. Since we are working to O (ci) in the 4t/ interactions, 
we will only keep one insertion of the £<*2 vertex in the functional potential operator (|58p . which thus reduces to 



y(4C,ci) 



6j ' <y t 



cim g 
2<?V 



d 3 zd 2 



diS 



s 



SJmn (5) 5jl M (z) 



diS 



SjQP iz) SjI q (z) 



(84) 



(where the notation (diS) /dj (z) implies that the partial derivative is acting only on the result of the associated 
functional derivative). The 2-point function is then obtained by taking derivatives of the generating functional (|51[) 
with respect to the matrix sources j and j\ i.e. 



(ul AB (x)U < , CD (y))= j- 



-5 



-5 



5jba (x) Sjl c (y) 



' S 5 ' 




8j' <5jt_ 


Z( 2U ) [0,0] 



(85) 



;t=0 



where the superscript (2U) indicates as in Eq. (|5"4")) that the corresponding quantity is evaluated by using only the 
2U contribution to the full soft-mode action. Eq. ([85)1 can then be rewritten as 



uIab (*) U<, C d (y)) - (ul AB (x) U^cd (y)) 



(2U) 



u 



t 

K.AB ' 



(a!) / tPzCW 14, (5) , Z7<, (5) 



U. 



<,CD 



(y) 



(2U) 



where 



Ul AB (a!) U <iCD (y) 



(2U) 



J, J' 



Sj BA (x)5jl c (y) Z(^)[0,0] 



SadSbcA 1 (y,x) 



(86) 



(87) 



7't=0 



is determined by the 2U part of the higher-derivative interactions only, while the 4U dynamics generates the pertur- 
bative O (ci) correction 



U 



<,AB 



(x) l <PzcW (z) u <lCD (y) 



(2U) 



C\m g 6 

W Sj BA (x)6 Jd . c (y) 



v [j,j t ] 
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which is due to the insertion of the vertex 

2<?V 



tl — 



cim g 
= 1 d 3 zd 2 



d 3 zC (4U) f — — 

<,C1 \5j'5jtJ Z( 2 ^)[0,0] 



^mjv (2) rfjt, (z) 



] 



(89) 
(90) 



(where d 2 acts only on the square bracket). 

After evaluating the functional derivatives in Eq. (|90[) . the contributions to the vertex v can be grouped according 
to the number of included source fields as 



z {2U) \j i f l 

v [ j,J ^ = z( 2 ^) [o, o] ( Urf + ^' t [ J ' ,J ] +w M t ^ t ]) 



(91) 



The constant Vd = v [0, 0] is the IR-divergent, disconnected part which originates from the "8" vacuum-bubble diagram, 
with all four legs of the vertex v [j, jt] closed among themselves, and reads 



Vd 



3, ,6 



-N a fj! 



ii(^c 1 )(27r) 3 S 3 (0) 



(92) 



(the integral 12 (£, c\) is defined in App. [B]). More generally, disconnected contributions to the n-point functions arise 
when all functional derivatives associated with external lines hit the Z^ 2U ^ Vd part of v The resulting 

products of free Green functions with Vd will play no role in our following discussion. All remaining terms require one 
(two) pair(s) (S/Sj) (S/Sjn to hit Vjp ( v jjUp) an d thus contribute exclusively to the connected 2-(4-)point function. 
The 2-line connected part of the vertex thus becomes 



%t [j, j f ] 



g 2 Nn 3 _ 



h (C,ci) 



d 3 : 



-jl N A- 1 d 2 A- 1 j NM + (fl(A- W) Am*' 1 d 



(93) 



(where the (?, ( d j) act on the first (second) argument of A 1 and integrals over the arguments of A 1 folded with a 
source are implied but not written explicitly). The 4- line connected part, finally, is 



[i,i f ] = / d 3 z(d l A~ 1 ]QN ) ] l Q A- 1 x 

x 2 (5, A-W) (jL/A- 1 ^,^) + (d 2 A- 1 j MP ) (jUA- 1 ^ 



(94) 



and contributes only to n > 4 point functions. (In both of the above expressions the "open" arguments (i.e. those 
not integrated over) are always the vertex coordinates z with di = d/dzi.) 

We now proceed with the calculation of the 2-point (and 2n-point) functions by evaluating the 4C7- vertex-induced 
functional 



n 



(2) 

ABCD 



-v [j,3^] 



3,3 



DC 



[0,0] 



■n 



(2) 

ABCD 



(95) 



which yields 



n 



(2) 



mh D [3,3^] (x, y) = 5adA c 1 b (y- x) (v^ + %t^t) 

+ J^"A„ <n a- (/ - e ) J*,*-* tt . ? )KD (?) fe + ,,,„„) 

*(%t [i,i f ] +^itjjt [i,i + ]) 



+ / dVA-My-^ico^) 
dV'i^ (/OA-^Z-f) 
a 2 (%-t [j, jt] +%tjjt [j, j f ]) 

^BA ^'Lc (?/) 



[i,j' f ] 



(96) 



16 



The 2-point function is obtained from II^ 2 ) by setting the sources j, to zero. Since then Vjj\ [0, 0] = v^\,^\ [0, 0] = 
and also the one-derivative terms as well as the term with two derivatives on Vjjjjp vanish, one is left with 

n% co [0,0](*,y)= / V ^:f\^ (97) 
ojba {x)6j' DC (y) 

(the right-hand side is independent of - hence one does not have to impose — explicitly). It thus remains 
to calculate 



J, J* , , 47V ff 4 [i 2 -d 
t = SadScb— —-5*2 & ci) d x j- 



SADScB-^-^h ci) d^A- 1 (y- x) (98) 



= -S AD S BC ^^A- 1 (x-y), (99) 

where we used the symmetry property (|76p and defined the "4[/-vertex-insertcd" soft-mode propagator 

x 1 o2 # A , 4oV _ , /" d 3 K K 2 6»(l-K 2 )e 4 ^ , 

as well as the dimensionless parameters 

Hence the c± correction (|88p to the connected 2-point function, which originates from the perturbative insertion of 
the 4U part (f53"| of the four-derivative interactions, becomes 

(uIab (2) / d 8 *^ (3 (y)) (2C/) = §^ n ABCD [0,0] (f,y) = -^fecglA- 1 (x - y) . (102) 

Our final result for the 2-point function is therefore 

(u<,ab (*) U <,cd (y)) = SadSbc A-^x-^ + ^A-^x-y) . (103) 

Its diagonal dependence on the group indices reflects the SU l(N c )x SXJr(N c ) symmetry of the soft- mode dynamics 
(HID. 



B. The 4-point function 



In the following section we sketch the analogous evaluation of the 4-point function. Although our main goal of 
computing c n corrections to the vacuum energy density involves (in our approximation) only the 2-point function, the 



explicit expression for the 4-point function 

■ (lll AB (zO U <iCD (z 2 ) rf <>EF (z 3 ) U <iGH (z 4 )) 

= ~ S ~ S ~ S ~ 6 y(4g,eQ 

Sjba (zi) 8^ DC (z 2 ) SjFE (z 3 ) 5ji H(J (Z 4 ) 

_. g(2U) _ 5 (AU,0{ Cl )) 

will be useful, too, because it allows 
density. Above, we have defined 

(2U) 



'5 5 ~ 


[j,jt] 


Sj' <Sjt_ 


Z( 2U ) [0,0] 



(104) 
(105) 
(106) 



useful, too, because it allows for consistency checks and an alternative estimate of the nonperturbative energy 
r. Above, we have defined 

t \ (: 
<,CD {z 2 ) Ul EF (z 3 ) U <fiH (z 4 )) 



[j, 

= SadScbA^ 1 (z2, zi) SehSgfA' 1 (z 4 , z 3 ) + SahSgbA^ 1 (z 4 , z{) SedScfA" 1 (z 2 , z 3 ) 



(107) 
(108) 
(109) 
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which shows that in the absence of AU interactions the 4-point function factorizes as 

[ U lab fa) u lcd (z 2 ) U[ EF (z 3 ) U LGH {z A )) = (u[ AB (zi) U LCD {z 2 )) (Ulef (30 U LGH (£4) 



+ ( u lab {zi)U lgh {za)) (uI ef (z 3 )U LC d (z 2 , 



and 



4 4t/ '° (ci)) = (uIab (*) U<,cd (z 2 ) J <PzC { Z\ (*) U* <tEF (z 3 ) U <tGH &)\ 



(2U) 



c x m q 



2g V Sj FE (z 3 ) 5j ] HG (z 4 ) Sj B A (zi) 5jl, c (z 2 ) 

ClTOg TT (4) ,-. -, _ _ , 

~ 2 ~ 2 ABCDEFGH ( Z 1,Z 2 ,Z 3 , Z4) . 



(110) 

(111) 
(112) 
(113) 



The above reordering of the functional derivatives reveals that Eq. (j9"5)l from the calculation of the 2-point function can 
be used as an intermediate step for the further evaluation. Indeed, to O (c\) the correction due to the AU interaction 
becomes 



n 



(4) 

ABCDEFGH 



(zi,z 2 ,z 3 ,za) 



SjFE (z 3 ) SjfjQ (z 4 ) 



n ABC d [i>i f ] (21,^2) 



3,P=0 



(114) 



which simplifies in the zero-distance limit of some of its arguments (as it occurs in our context) and then renders, 
together with Eq. f|109[) . the coordinate dependence of the 4-point function explicit. 



VI. SOLUTION OF THE GAP EQUATION AND PHASE DIAGRAM 

As explained in Sec. IIVB| the above expressions for the 2- and 4-point functions are to be evaluated at the saddle 
point S =: (/x£) of the integral over the auxiliary S field, i.e. at the minimal- action solution of the gap equation 

(x)U <t Bc(x))=S AC . (115) 
After inserting the result (|103[) for the connected 2-point function (and recalling Saa — N c ), Eq. (|115l) turns into 



U <,AU 



= 1. 



(116) 



The zero-distance limits of the propagator ([73")) , (|74l) and of the AU correction (|100p are then expressed in terms of 
the integrals i n , j n (defined and evaluated in App. |Bj as 



A- 1 (0) = ^/i 3 ? 2 fec 1 )j 2 (£ >Cl ) 



(both are of O (g 2 /C))j so that the gap equation assumes its final form 



«i (€,ci) + 2ct-h (C,Ci)j 2 (£,ci) 



= 1. 



(117) 
(118) 

(119) 



(For ci = and £ = 1, Eq. (|119|) reduces as expected 17] to 711 (£) = 1 (except when 7 diverges, see below).) 
This equation and the £ dependence of the integrals (|B1|) . (|B2[) render the nonperturbative character of the solutions 
£ (/i, c±, C) explicit. As already mentioned, it reflects the infinite subset of diagrams required to generate a finite mass 
gap. 

Before finding the solutions of Eq. (II as functions of the variational parameters, we specialize £ to the form 
required by continuity of G~ l (k) at k = fi (ci. Eq. (f2U|) . 



(120) 
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and adopt the one-loop Yang-Mills coupling 



7 0) 



9YM N C iVe=3 



24 



11 In 75^- 



(121) 



to render the \i dependence of 7 explicit. The coupling QymSm) * s i us ^ H/10 times larger [32|,[35| than the coupling ([57)1 
obtained from integrating out the high- momentum modes [351 ] governed by the UV covariance (jlip . (This discrepancy 
can be mended by admitting non-transverse modes to the UV covariance [II].) We prefer to use the Yang-Mills 
coupling because it facilitates direct comparison with the numerical results of Ref. 17| . 

The above considerations imply that the solutions of Eq. (|119j) depend on just two variational parameters, i.e. 
the RG scale /i > which enters through the coupling (|121[) and ci < 1 which controls the leading momentum 
dependence of the infrared covariance G^ 1 (k). Our next task will be to determine the critical line /i c (ci) in this 
parameter space, i.e. the subspace which joins the points where the phase transition takes place and where thus the 
(dis-)ordcr parameter vanishes, 



£(Mc(ci),ci) = 0. 



(122) 



This line can be found analytically since for £ = the integrals i\ , 12 and j'2 in the gap equation (|119p simplify to (cf. 
A PP . E2J 



1(d) := i\ (0,ci) 



arctanh «/ci 



and 



h (0, Cl ) = i 



1 



h (0, ci) 



*(<*)-! 



C\ 



1-Cl 

In order to exhibit the [i dependence, we first resolve the £ = gap equation for 7 (/1), which yields 

1 

72" 

Set 



+ *(t-i) + — 7-1 = 0. 

1 - Cl / Cct 



(123) 



(124) 



(125) 



The two solutions of this equation can then be resolved for /i and combine into the critical line /i c (c%). With 7 (/x) 
from Eq. (|12ip it takes the explicit form 



Mc,l,2 (£1) 

Aym 



exp 



48 



(l-ci)[l-*(ci)] [l + (l-d)*(ci)] 



11 



(1 - ci)*( Cl ) ± ^5? (ci) (1 - Cl ) 2 _ 4(1 - ci) [1 - ci?(ci)] 



(126) 



In order to facilitate the discussion of the vacuum phase structure encoded in Eq. (I126[) , we plot the phase boundary 
in Fig. [2] It reveals that Eq. (|122[) is satisfied by two values of /j, c for each c\ for which a solution exists, except at 
the maximal and minimal values where \x c becomes unique. Vice versa, there are two ci for each fi c for which Eq. 
(I122p holds. The critical line (|126[) covers the limited parameter ranges 



0.5 < 



A 



< 



YM 



and 



0.48 < ci < 1. 



(127) 



(128) 



This prevents the minimal- energy solution £* (to be determined in Sec. IVIII Al) from attaining unacceptably large 
values of \x and |ci|. The normalizability condition c\ < 1 is automatically satisfied in the existence region of gap- 
equation solutions. Equation (|126p further shows that nontrivial (i.e. nonzero) solutions of the gap equation exist 
only when the gauge coupling exceeds a critical value, i.e. for 



9 2 (M) > 9c (ci) , 

as expected on physical grounds. The maximal critical coupling corresponds to ci = 0: 

^2 



at (ci) < gi (ci = 0) 



(129) 



(130) 
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Only part of the critical line (|126[) is physically trustworthy, however. Indeed, its validity range is limited to |ci| <C 1 
for which the O (c\) evaluation of the AU correction included in Eq. (|126[) is accurate, and to /x/Aym ^ 5 — 6 which 
justifies the O (g) approximation. These two restrictions eliminate most of the /z Cj i branch and the lower-/z part 
of the /i Cj 2 branch. (Although the c\ = solution with /i c /Aym —> 1 and thus g 2 (/i c ) — > oo invalidates the 0(g) 
treatment, it nevertheless demonstrates that the gap equation does not reduce to 711 (£) = 1 for c\ — > if 7 (p c ) 
diverges simultaneously.) 

In Fig. [21 finally, we plot the numerically generated solution £ (/j,, c%) of the gap equation (II 10[) in the physically 
reliable parameter range. As expected, the £ > region is surrounded by the critical line (I126[) and exists for 
0-5 i$ /VAym S 8.86. In addition, Fig. [3]reveals that the (dis-)order parameter goes to zero continuously, i.e. that 
the transition from the disordered to the ordered phase is of second order [6{| . 



VII. VACUUM ENERGY DENSITY 



We have now assembled all the infrared-mode information necessary to evaluate the expectation value (f2"5j) of the 
Yang-Mills Hamiltonian density 

Hym=1(E?E? + B?B?) (131) 

in the generalized trial vacuum state ([T]). (We have adopted the Aq = gauge, and defined the chromoclcctric 
and -mag netic fields as Ef = F& and Bf = \e ljk Ff k .) Of course, the Hamiltonian (| 13 1[) is formal and needs 
renormalization. Since we are interested in its matrix elements between Poincare-invariant vacuum trial states with 
perturbative one-loop corrections only, a regularization of the vacuum matrix elements by a momentum cutoff Auv 
will be sufficient [TtJ • In fact, the Auv dependence of the vacuum energy density (l2"5t can be completely removed by 
normal-ordering the Hamiltonian til 31[) with respect to the noninteracting vacuum. This amounts to subtracting the 
ultraviolet-divergent energy density 



(Hym)o = 2 (N? 1) = 2 (N? - 1) J ^6 (A^ v - P) ^ = %^A UV (132) 

(where we wrote — k for the free gluon energy and reinstated ft) of the noninteracting vacuum, i.e. the sum over 
the zero-point energies of two transverse, massless vector modes, from (Hym)- (Recall that the momentum cutoffs are 
imposed on integrals over the gauge-invariant U modes and therefore do not compromise (residual) gauge invariance.) 



A. Hard-mode contribution 



According to our strategy for evaluating the trial energy density (Hym) outlined in Sec. IIII A[ the first step consists 
of integrating over the static gauge fields and the hard modes U > . In the following section, this will be done for the 
chromoelectric (-E 2 ) and -magnetic (-B 2 ) contributions separately. We start from the intermediate matrix element 



{{(Ef{x)E]{y)))) 





>" 


iS i8 
SAf(x) SA^y 


yV'o 


'A 


j DAtp 


A u 


V- 


A 





S^G- 1 (x,y) 



d 3 Zl / d 3 z 2 G- 1 (x-z 1 )G- 1 (y-z 2 )(((A-(z 1 )A b J (z 2 )))) 



(133) 



(Ef — i5/8A®, ci. App. [5J for the chromoelectic 2-point function. With the help of the expression (|A21j) for the 
gauge-field matrix element in terms of the functionals a [U] and M. [U] as defined in Eq. (|A14j) . this specializes to 



(((E 2 ))) := (((Ef (x) E? (x)))) = 3 {N 2 l) G~ l (x, x) 

- [ d 3 z 1 d 3 z 2 G- 1 (x-z^G- 1 (x-z 2 ) [M^ laa (z u z 2 )+at (z x ) a? (z 2 )] 



(134) 
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The analogous expression for the chromomagnetic contribution involves 2-, 3- and 4-point functions of the gauge field. 
After rewriting them with Bf = e ijk (d 5 A% + § / obc A^) and Eqs. ((Ml"]) - (|X23| . one finds 

{{{B 2 ))) := ({{Bt (x) Bf (a?))» - £^^ m d^d^M^T C?!^)^ + W ^-a^< 
+ [X- 1 (af, xt n d 3 a a k + d zu ,M- x {A^TL of 



(135) 



(where all spacial arguments are set to x after taking the derivatives). 

Now we compute the hard-mode contributions to the chromoelectric and -magnetic vacuum energies, as outlined 
in Sec. IIII Al and App. [5] This amounts to evaluating the intermediate matrix elements (|3"Tj) for E 2 and B 2 , i.e. 



((E 2 ,B 2 )) 



J Dcj) [(((E 2 , B 2 ))) + Q (g)] exp {-r b ,> [0]} 
/D0exp{-r b ,> [</»]} 



(136) 



to leading order in g. The hard-mode action r b .> [0] is defined in Eq. (|A42[) . Since it is bilinear in </>, the integrals 
in Eq. (|136[) are Gaussian and can be calculated analytically, cf. Appendix [3] This yields 



«£ 2 » = (^c 

1 



d 3 k 

(2^ 



' (A uv - k 2 ) G- 1 (k) 



N 



d 3 k 



2 )G- 1 (fc) + «^ 2 )) u< +0(. 9 ) 



which includes the contribution 

(( E2 ))u < -~Jg2 

from the soft modes U < , and 



d 3 z, / d 3 z 2 G~ 1 [x-zt) G- 1 {x-z 2 ) {{L a <ti {zx)L a <ti (z 2 )» 



{{B 2 )) = (N 2 1) 



d 3 k 
(2^ 



9(A 2 JY ~k 2 ) k 2 G(k) + 0(g). 



(137) 



(138) 



(139) 



With 9 (k 2 - n 2 ) G(k) =9 (k 2 - ^ 2 ) fc" 1 (cf. Eq. (HU) one finds the Auv dependence of {{E 2 )) and {{B 2 )) to be 
identical. Hence the lowest-dimensional gluon condensate (-F 2 ) = 2 [(-B 2 ) — (-E 2 )] (to be discussed in Sec. IVIII B[) 
is Auv-hrdependent by itself while the Auv dependence of ((Hym)) has to be removed by renormalization. As 
noted above, at our level of approximation this amounts to normal-ordering the Hamiltonian, i.e. to subtracting the 
non-interacting vacuum expectation values (-E 2 ) = {B 2 ) Q = (N 2 — l) G m L (x, x). The results are 



((:E 2 :)) = (N 2 -l) 



d 3 k 



and 



((:B 2 :)) = (N 2 c -l) 



:9U 2 



d 3 k 

(2^ 



G" 1 (fc) - k 



((E 2 )) v< +0{g) 



(v 2 - k 2 ) [k 2 G (k) -k]+0 (g) . 



(140) 



(141) 



These expressions are indeed Auv-independent and exclusively receive contributions from momenta k < /.i. (Both 
properties have to become manifest at this stage because they are unaffected by the remaining integration over U < .) 
Moreover, the above results provide a nontrivial check of the fact that by restoring residual gauge invariance the 
integration over </> has removed the longitudinal hard-mode contributions. 



B. Soft-mode contribution in the disordered phase 



It remains to calculate the chromoelectric soft-mode contribution (|138[) to the vacuum energy density. Here the 
impact of the generalized IR covariance (1131) becomes fully nonperturbative and the dependence on the variational 
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parameters enters partly through the solutions of the gap equation. In order to evaluate Eq. (jl38[) . we first perform 
the remaining part of the U < integration (cf. Eq. (|33l) ). i.e. we integrate over £, which yields 

(E 2 ) v< = ~ J d 3 z, J d^z.G- 1 (x- z,) G- 1 (x-z 2 ) (L%. (%) L% ti (z 2 )) . (142) 

(Recall that in the mean-field approximation of Sec. IIVBI the £ integration amounts to evaluating 
((£< } i (zi) i< j {z 2 ))) at the saddle-point solution £ (//, ci) = (/i, ci)] of the gap equation (|119l) .) We further note 
that the integrals over G _1 (x — z\) G _1 (x — z 2 ) in Eq. (11421) have most of their support at distances \zi — z 2 \ < 
where U < (z 2 ) {/< (z\) ~ 1. This allows us to approximate 



(L% ti (fi) L% ti (z 2 )) = 2 (tr {diUi (z 2 ) U< (z 2 ) u\ (*) d l U < {£,))) 

2(tr{d l uUz 2 )d l U < {z 1 )}) 
(as in the 2U approximation (|4"4")l to the Lagrangian (|43l) ). From Eq. (I103P one then has 



(L« i (z 1 )L a < Jz 2 ))~2N% 



A" 1 (z 2 - A) + ^A" 1 (z 2 - z x ) 



C/1 3 



(143) 
(144) 

(145) 



and with Eqs. ([T3]l. ((74)) (specialized to c„> 2 = 0) and (11001) further 



(L'sWL'j (z 2 )) ~ 47V c7 r 2 ^ 2 



rf 3 « 

(27T) 3 



k 2 (1 - k 2 ) e «M*>-*0 
1 | 2C17C-H2 (C.ci)k 2 

K 2 (l- Cl K 2 )+£ 2 [ K 2(l- ClK 2) + ^2] 2 



which has the coincidence limit 



(L",, (*)£",< (z))^27V c ^ 2 



h (£,c x ) +2ci-? 2 (C,ci) j 3 (£,,ci) 



(146) 



(147) 



With the explicit expression (p~3|) for the generalized IR covariance (with c;> 2 — 0) Eq. (jl42[) evaluates further to 



1 -Q\ 



1 P 
1 - Ci — 



X / d 3 0! 

and, after inserting Eq. (|147|) . assumes its final form 



cfz^Z-^e^-^ (L a Ki (zx) L a Ki (2b)) 



2tt 2 



7' 



„ 7 /- o~ 

«2 - 2cii 3 + cp 4 + 2ci-i 2 [J3 ~ 2c x j 4 + cjj 5 



(148) 



(149) 



where the integrals i n , j n are understood to be evaluated as i n Ci) ,ci) , j„ (£ (m? c i) : c i)- (Note that these 
integrals receive most of their contributions from vacuum modes with momenta k *~ [X (cf. App. lAj).) 

The rather complex parameter dependence of Eq. (|149l) simplifies considerably for small |ci| <C 1. Indeed, after 
specializing the IR mass parameter £ to £ct (ci) = (1 — ci) = 1 + ci + O (c 2 ) (which ensures a continuous covariance, 
cf. Eq. (|120p 'l and using the small-ci expansions (|B66[) . (|B67I) of the integrals i n and j n , the soft-mode contribution 
eu < — (E 2 ) u 1 2 to the energy density becomes 



iV 2 , , 

c 4 f ~ 



4tt 2 



{?2 ® + CI [? 2 (I) - 2? 3 (I) + J4 (I) + 2 7 ? 2 (I) j 3 (I)] } + O (c?) 



AT 2 

-^ 4 (1 



(150) 
(151) 



for |ci| -C 1. (The last expression provides a very reasonable approximation for |ci| < 0.1.) With ci = and m g = fi 
(and approximating TV 2 — 1 ~ N 2 ), finally, the above expression reduces to the result 



N 2 N 2 
e v< ( M , Cl = 0, C = 1; = -^h ® = 



1 



1 



- - £ ( 1 - £ arctan ^ 



(152) 



of Ref. [17J, as it should. 
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C. Energy density in the disordered phase (i.e. for ji <C Auv) 



The (normalized) integral over the auxiliary field S turns the matrix element ((Hym)) into the trial vacuum energy 
density (Hym) (cf. Eq. (I3"3")) '). This integral is nontrivial only for the U < and hence £ dependent part (|138p of the 
integrand (JT3TJ) . Its evaluation in the previous section yielded Eq. (|149p . Combining these results and separating the 
complete vacuum energy density e = E/V = (Hym) into hard and soft contributions, 



(/", ci, C; i) = (Hym) = £> (m) + £< (/x, ci , C; f) 



(with £ still unspecified), one finds 



1 



£> (M) = 2 (A 2 - 1) 



and 



£< (M: c 1)C;I) 



AT 2 



rf 3 fc 

(27T) 3 
d 3 fc 



(A^ v - fc 2 ) - 9 ( M 2 - fc 2 )] [G- 1 (k) + fc 2 G> (fc)] 



(27T) 



(m 2 - fc 2 ) 



G" 1 (fc) +fc 2 G< (fc) 



+ o( i5 V (^ci,C;e) 



The energy density e> of the hard modes involves only G> (k) = k = k 2 G > (fc) and evaluates further to 

[9 (Al v - fc 2 ) - 6 („ 2 - fc 2 )] k - = ^ (A uv - m 4 ) , 



£> ( M ) = 2 (7V C 2 - 1) 



(27T)' 



(153) 
(154) 

(155) 
(156) 



which is the (regularized) zero-point energy density of the two transverse, massless vector modes with energy uj (fc) = fc 
(recall that the integration over <f> has removed the longitudinal-mode contribution) and reflects the built-in asymptotic 
freedom of the fc > [i modes. As anticipated, the subtraction of the free vacuum energy density (I132p in the course of 
normal-ordering cancels its Auv dependence. 

The soft-mode contribution is mainly of nonperturbative origin and therefore structurally more involved. Inserting 
the covariance (|20|) and its inverse (f21"]l (for c n >2 = 0, as before, and C = m gl v) into Eq. (|155|) yields 



4tt 2 



c 



3ci 



1 



c\ arctanh y/ci 



(157) 



so that the total energy density (|153l) . after dropping the constant zero-point contribution (A^ 2 — l) A uv / (87r 2 ), 
becomes 



1M 4 



4ir 2 (ci 



3 (3Cci - 5C + 5) Ccf + lOci + 30 arctanh ^/cT 



30 



(158) 



After further approximating A 2 — 1 ~ A' 2 , ensuring the continuity of G< (fc) by imposing Eq. (I120p and using Eq. 
(11491) . we obtain our final expression for the energy density in the disordered phase, 

e (fj,, ci) := e (p, Cx,(ct (ci) ; £ 0", ci)) 

_ N 2 ^ 4 [4c? + 10c 2 - 50ci + 30 1 - c x arctanh jc{ 



4tt 2 ' 



30cf(l-ci) 



i 2 - 2ci« 3 + cfu + 2jd (1 - ci) h (ja - 2cij 4 + c 2 j 5 ) 
1^ 



(159) 



where the integrals i n , j n are evaluated at £ (/x, ci). For c\ — > 1 several terms of the energy density (|159[) diverge and 
prevent the vacuum from encountering the limiting instability. For |ci| <C 1, on the other hand, the rather complex 
Ci dependence of the full vacuum energy density (|159[) simplifies to 



e(p, c x ) 



N 2 A 



Aw 2 



+ u (£) - 2? 3 (£) + 2 7 ? 2 (£~) j 3 (£~) 



cx\ + (cf) 



(160) 



to the energy density of Ref. 



(where the solution £ of the gap equation is evaluated at ci = 0). For ci = and £ = 1, finally, it reduces as expected 

• m 

N 2 

e (n, cx = 0, C = 1; £ 0, cx = 0)) = 



h £^ [ 1 - £ arctan i 



15 



(161) 



The full energy density (II 50[) is plotted in Fig. 2] and will be discussed further in Sec. IVIH Al 
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D. Energy density in the ordered phase (i.e. for /i ;§> Aym) 

As argued above, the mean-field approximation is reliable in the disordered phase, i.e. for those /i < [i c < 8.86Aym 
(cf. Eq. I127[) for which £ is not too small. In order to get a more complete picture of the vacuum energy density 
and its \i dependence, however, one has to evaluate the soft-mode contributions in the ordered phase as well, i.e. for 
/i 3S> Aym where the mean-field approximation breaks down. However, in this phase the Yang-Mills coupling becomes 
small [70] . i.e. 

g 2 fa ^> A YM ) < 1, (162) 

so that the soft-mode dynamics can instead be treated perturbatively. After restricting ourselves in accordance with 
our previous approximation scheme to the 2U contributions in Eqs. (|42|) and (|44|) (again with Cj>2 = 0), this dynamics 
becomes 



(163) 



Since in the ordered phase fluctuations y?" around U < ~ 1 are small, one can furthermore truncate the weak-coupling 
expansion of U <1 i.e. 

U< = exp (w<A a ) = 1 + ig<A a + O (g 2 ) , (164) 
which to O (g 2 ) yields (with tr {A a A fc } = 2S ab and denoting the Fourier transform of ip as (p) the bilinear action 

1 f d 3 k 



(165) 



From this "kinetic" term one reads off the static k < fj, propagator of the <p modes (which contains, in contrast to the 
soft-mode propagator ([75]). ([74]) in the disordered phase and as a consequence of Goldstone's theorem, no mass term) 
as 

xab r j3jL -ife(zi-z 2 ) 

{v% ma m _ i_ J At 9 p _ e) C W - WF . dee) 

The correlator (|143| . which determines the chromoelectric vacuum energy density due to the soft modes, then becomes 
(with S aa = N 2 - 1) 

(L% ti fr) L a <ti (z 2 )) ~ 2 (tr [d^ (*) d i U < (z 2 )}) ~ Ag 2 (d^ (z\) (z 2 )) (167) 

= 2 (A 2 - 1) tl. [ ii^e-W-^). (168) 
V ; ( i (2tt) 3 1-ci^ 2 

After inserting this result into Eq. (|148j) one obtains 

= - (^c - 1) ^ / |^ (1 - « 2 ) (1 - rW) (169) 
A 2 - 1 ,, fl 1 



4tt 2 ,r V3 5 
Hence Eq. (| 157[) yields the infrared contributions 



Cm 4 U " 7 c i • ( 17 °) 



A 2 - 1 zz 4 f~3C 2 4 + 5C 2 c? - 5ci - 15 arctanh v^r". 



to the energy density which specialize with Cct (ci) = (1 — ci) 1 to 
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Together with the hard-mode contribution (|156|) (which is identical in both vacuum phases) and after discarding the 
zero-point contribution, finally, the total vacuum energy density in the ordered phase becomes 
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(174) 



Several properties of Eq. (|173[) are noteworthy. First of all, in the C\ — > limit it reduces, as expected, to the result 
ofRef. [nl (cf. their Eq. (4.24)), 



£ (n,Ci = 0) 
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(175) 



The total energy density diverges for c\ — > 1 where the vacuum wave functionals would become unnormalizable, 
furthermore, which prevents c\ from growing beyond unity during energy minimization. Most importantly, however, 
for /i ^> Aym the negative hard-mode contribution to the energy density (i.e. Eq. ()156[) without the Ay V term) 
is overcome by the positive contribution (|172[) from the soft modes. Indeed, the (large /i) vacuum energy density 
(11731) is positive for all c\ < 1 and hence a monotonically increasing function of fi. In addition, Eq. (jl 73(1 increases 
monotonically with ci in the range — 2 < ci < 1 which includes those |ci| <C 1 values for which our perturbative 
O (ci) treatment is reliable (and thus our variational result c\ ~ 0.15 to be determined in Sec. IVIII A|) . Under these 
circumstances it is reasonable to expect that the perturbative result (| 1 73[) remains at least qualitatively reliable for /i 
values down to the phase transition [l7| at fi c < 8.86Aym- 



VIII. VARIATIONAL ANALYSIS 



We are now prepared to minimize the vacuum energy density in our trial-functional family according to the Rayleigh- 
Ritz variational principle. Above we found that (inside the validity range of our approximations) the energy density 
(I159p in the strongly-coupled disordered phase decreases monotonically with increasing fj, up to the phase transition, 
whereas its counterpart (|173[) in the ordered weak-coupling phase monotonically increases both with fj, and c%. The 
combination of these results indicates that the vacuum energy density attains its minimum at the phase boundary 
in the disordered phase, where the number of massless particles becomes maximal. Indeed, the number A^ 2 — 1 of 
massless Goldstone modes in the ordered phase (at very large \x all other modes have m ~ Auv and are too massive 
to contribute significantly) doubles at the transition where they are joined by degenerate parity partners. Hence the 
massless particles generate roughly twice the internal energy at the second-order transition, which becomes maximal 
at the critical /i* [13] • (This argument is not affected by the higher-gradient interactions (|43|) since those share the 
symmetries of the leading term (02]).) 

Our program for the following subsections will be as follows: after determining the energy minimum at the boundary 
of the disordered phase quantitatively (Sec. IVIII Aj) . we evaluate the associated, four-dimensional gluon condensate in 
Sec. IVIII Bl The resulting vacuum field distribution and its physical interpretation are discussed in Sec. IVIII Cl and 
its impact on the phase structure of the soft- mode dynamics is subject of Sec. IVIII Dl In Sec. IVIII El we comment on 
the qualitative role of the higher-gradient interactions in the soft-mode Lagrangian (|43[) . Finally, in Sec. IVIII F I we 
mention several options for improvements and future applications of the gauge-invariant variational framework. 



A. Vacuum energy minimization 

In Fig. @]the energy density (|159l) of the disordered vacuum phase (calculated on the basis of the numerical solution 
of the gap equation, cf. Fig. [3]) is plotted as a function of [i and c\. The plot range of fi includes those regions in 
which (i) the one-loop evaluation of the hard-mode contributions remains reasonably accurate (corresponding roughly 
to fj, > 4Aym), (h) the system stays in the disordered phase where the nontrivial solution £ of the gap equation exists, 
and (iii) the mean-field approximation remains approximately valid. Together with Eq. (|127p these conditions require 

4A YM < M < 8.86A YM - (176) 

For illustrative purposes the plot range of c x is chosen to cover almost the full existence range —0.48 < c x < 0.95 
of gap-equation solutions, on the other hand, although our perturbative treatment of the 4U interactions becomes 
questionable close to the upper limit. 
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The arguably most prominent qualitative feature of the vacuum energy density e (/x, c\) in the disordered phase is its 
monotonic decrease with both /i and C\ which continues until the nontrivial solution of the gap equation ceases to exist 
at the second-order phase transition (cf. Fig. HJ. This essential feature manifests itself already in the linearization 
(|160l) of the Ci dependence around c\ = 0. As a consequence, the energy of the disordered vacuum is minimized 
at £ = 0+, i.e. at the disorder-order phase transition, for each admissible value of c\. In order to find the precise 
minimum of e and the corresponding parameter values fi* and c*, we plot the energy density along the critical line, 
i.e. 

e(c x ) := e(fi c (c x ),ci) , (177) 

in the range ci € [—0.48, 1] where £ (/i c (ci) , ci) = 0. As can be read off from Fig. [51 the minimum of e (ci) is attained 
at 

c{ ~ 0.15 with fi* = fi c (c{) = 8.61A YM - (178) 

Evaluating qym (m*) = 9ym (/-**) / (^ 7r ) — 0-27 (to one loop) at this scale confirms that the running coupling remains 
small enough to justify the perturbative treatment of the hard modes. (It exceeds ceyM (Mc (ci = 0)) ~ 0.26 by less 
than 2%.) Similarly, the resulting cj < 1 justifies our perturbative treatment of the 4U contributions. In addition, 
the rather small curvature of the critical line (|126[) around its maximum 

/4 roax) = Mc (ci = 0) = cxp (^j ~ 8.86A YM (179) 

at ci = causes our dynamical IR mass scale \i* ~ 8.6Aym to be only about 3% smaller than /ic inax ' ) . This has 
important consequences for the value of the gluon condensate to be discussed in Sec. I VIII B I and for other vacuum 
scales. The corresponding IR gluon mass, e.g., turns out to be 

(C*) 

1 



m* = ~ 10.14A YM . (180) 



(For Aym — 0.15 GeV [13] this value is somewhat larger than the potentially related mass ttilg ~ 1.1 GeV extracted 
from the intermediate-momentum behavior of the Landau-gauge gluon propagator on the lattice (second reference of 
34]).) 

The quantitative improvement of the vacuum description due to our generalized trial functional basis can be 
measured in terms of the achieved variational bound on the Yang-Mills vacuum energy density. Comparing the 
minimum of our trial energy density (|159[) , 



e(c*) ~ -210.59A YM , (181) 

to the value E{c\ — 0) ~ — 187.52A YM generated by the lowest-order covariance G^ 1 (k) = [i shows that the Ci 
corrections reduce the vacuum energy density by about 11%. Hence variational optimization in our extended trial 
space has produced a rather substantial improvement of the vacuum functional. (The massive gluon propagator (|24l) . 
on the other hand, is energetically disfavored: its IR behavior is approximately reproduced by ci = —1/2 and yields 
e(c\ = —1/2) ~ 0. We expect sizeable corrections to the O(ci) treatment of the 4U interactions for ci = —1/2, 
though, and further note that Eq. (|2"4")1 does not match to the UV covariance (| 1 1 1) at k = fi, in contrast to our IR 
covariance.) 



B. Gluon condensate 



Gluon condensates, i.e. the vacuum expectation values of gauge-invariant local operators composed of gluon fields, 
are among the key amplitudes which characterize the Yang-Mills ground state. The most important gluon condensate, 
which dominates the power corrections in the operator product expansion e.g. of the glueball correlators (46| . is the 
expectation value of the lowest-, i.e. four-dimensional gluonic operator F 2 = F^ l/ F a ' fiL ' — —2 (EfEf — BfBf). In our 
trial state, and to lowest order in the gauge coupling, this condensate becomes 



(F 2 ) = 2 [(B 2 ) - (E 2 )} 
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where we used Eqs. fjl3T[) - (|139|) and again performed the X integrals in the saddle-point approximation. As 
anticipated above, the hard-mode (i.e. k > /j) contributions to the chromo-electric and -magnetic parts cancel, 
which indicates that the gluon condensate (|183|) is renormalized at ix. After inserting the IR covariance (|20|) and the 
propagator (|2"Tj) (with c„>2 = 0, as before) the condensate becomes 
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(184) 
(185) 



At the border of the disordered phase, where £(a< c (ci) , ci) — and where the energy becomes minimal (cf. Sec. 
IVIII Ap . we furthermore have from the £ — > limit of Eq. (|149l) and the expressions for the integrals i n (0,ci) and 
j n (0, ci) given in App. IB 21 that 
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Combining the above results and specializing as before to £ c t = (1 — c\) then results in our final expression 
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(7* = g 2 (/i*) N c /ir 2 ~ 1.012) for the gluon condensate. For small ci Eq. (|187p expands into powers of ci as 



<F 2 ) = 



N 2 -l 
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(186) 



(187) 



(188) 



which shows that (-F 2 ) grows with Ci when ci 1. In the uncorrected case, with ci = and the corresponding value 
a Y M (a*c (ci = 0)) /7r = 1/ (4/V c ) or 7* (/i c (c c = 0)) = 1 (cf. Eq. (jlOip ) for the gauge coupling as well as Eq. (|179p for 
the IR scale, and furthermore setting N 2 — 1 ~ A 2 , N c = 3, one finds [13] 
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Our improved value of the gluon condensate, on the other hand, is obtained by inserting the energy minimizing values 
Ci ~ 0.151 and = 8.606Aym from the last section and the corresponding coupling 7* ~ 1.012. The result 
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2X = 20.87A^ M 



0.011 GeV 4 



(190) 



is about 25% larger than the uncorrected value (|189[) . In the second equation above we have specialized to Aym — 0.15 
GeV, which allows for comparison with the value range ((a/71") F 2 ) = 0.0080 — 0.024 GeV 4 obtained from QCD sum 
rules (for references see e.g. Ref. [Igj]). (We refrain from using the currently preferred value for Aqcd which contains 
quark contributions.) Our condensate value prediction (I190[) lies comfortably within this standard range. 



C. Vacuum field distribution and quasigluon velocity 



Our variationally optimized wave functional determines the distribution of the gauge-field modes Ai (k) and thus 
contains new information about the field composition in the vacuum state. In our approximation, the IR covariance 
(|20| with c„> 2 = and Eq. (J23J) takes the form 



G- 1 (k) 
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1-ci 



1-ci- 
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(191) 



By construction, the covariance (I19ip is non-negative in the physical parameter range fj, > 0, c\ < 1 and thus yields a 
normalizable wave functional. In addition, G^ 1 (k) becomes larger (smaller) for positive (negative) values of c±. This 
holds for each k < /1 independently, but gets more pronounced for smaller k (cf. Fig. [T]). Hence the Gaussian weight 
factors 



exp 
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which (when multiplied for all k < /j,) make up the IR part (|10l) of the unprojected core functional, become narrower 
(wider) in Fourier-mode space [7l| . As a consequence, large Af (k) contributions to the IR part (flU)) of the unnormal- 
ized Gaussian vacuum functional, and thus to the functional integrands of any amplitude, are increasingly suppressed 
(enhanced) towards smaller k £ [0,/x]. 

Hence our positive result c* = 0.15 reshapes the vacuum field population in a specific way: relative to the uncorrected 
case ci = the attractive IR interactions generate a depletion of the ultralong-wavelength k — > modes and an 
enhancement of the k ~ fi modes. This effect helps to prevent the Savvidy instability which constant chromomagnetic 
fields experience in the Yang-Mills vacuum [47} . It also contributes to generating the expected average wavelength 

A ~ Ayl (193) 

of the vacuum fields. In contrast, the massive vector covariance (|24|) with c\ = —1/2 would populate the medium-soft 
modes with k < /j, less strongly (cf. Fig. [1} although this is energetically disfavored. (When comparing the overall 
weight of the k < \x modes associated with different covariances, the G~ 1 dependent normalization factor J\Tg has to 
be taken into account as well.) Since the higher-derivative soft- mode interactions in Eq. (l43l) do not affect the k = 
modes, the above findings underline the importance of finite-momentum IR modes in shaping the vacuum structure. 

Regarding Eq. (|191[) as the dispersion relation w (k) of "quasigluon" modes in the vacuum, one may further relate 
the parameter c\ to the dimensionless quasigluon group velocity 



v(kj 
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Since our approximate dispersion relation (|191j) is quadratic in the momentum, one may further define a (momentum- 
independent) effective kinetic gluon mass 

- - 1 C V (195) 



2cim g 2ci 

which relates velocity and momentum as k — m g v. Note that m g is negative for > c% > 1. With v (k) := v (jtj 
one further has 

\c x \ = ^-v{n) (196) 
2m g 

which shows that c\ is related to the quasigluon group velocity at the IR renormalization point k — \i. Specializing 
to a continuous covariance by identifying fi/m g = 1 — c\ then yields 

This relation implies among other things that the velocity diverges toward the bound c\ — > 1, and that the admissible 
velocity range is restricted to oo > v (/z) = 2 |ci| / (1 — c\) > —2. Our result c* = 0.15 corresponds to v* := v (fj,*) = 
0.35. 

The positive sign of c\ implies that the effective kinetic mass (|195[) is negative, i.e. that the quasigluon velocity in 
the vacuum is opposite to its momentum. Hence these quasigluons decelerate when an external force is applied, in 
stark contrast to the behavior of free gluons. (An interesting and related issue is the vacuum response to external 
electric and magnetic background fields, cf. e.g. Ref. [27}.) In this respect our IR dispersion relation is reminiscent 
of optical phonon branches in periodic structures over a Brillouin zone, or of the dispersion of electrons moving in 
such structures. In the latter example, the effect is generated by the lattice exerting a large retarding force on the 
electron that counteracts and overcomes the applied force. Of course, in our aperiodic momentum space no Bloch- 
type oscillations are generated by external color-electric fields. Nevertheless, quasigluons (with their small scattering 
amplitudes) or quarks could show a negative differential resistance, in remarkable contrast to the impact of a simple 
infrared mass term of the type ((24)) . The origin of these effects in the Yang-Mills dynamics should be investigated 
further. On a cautionary note, however, one has to keep in mind that not all properties of the gauge- variant core 
functional © play a physical role (cf. Sec. Ill Bl and Ref. [20}l. 



D. Phase structure and transition order 



It is instructive to analyze some additional features of the c\ ^ interactions which have direct impact both on 
the phase structure of the soft-mode dynamics (|42[) . (|43[) and on the vacuum energy minimum. To begin with, we 
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recall that the vacuum instability for c\ > 1 is encoded in c\ — > 1 singularities of the energy density and in zeros of 
the amplitude entering the gap equation. Hence the constraint C\ < 1 does not have to be imposed by hand but is 
automatically enforced by the diverging vacuum energy density and by the nontrivial saddle-point solution ceasing 
to exist. Since the energy density (|159[) decreases monotonically with /i, furthermore, its variational minimization 
could drive /i far below the validity range of the perturbative integration over the hard modes (cf. Sec. IVII A [I . were 
it not stopped at a sufficiently large fi* by the vanishing of the gap-equation solution ci), i.e. by the order- 
disorder phase transition. This behavior generalizes to the c\ dependence. Indeed, although the energy density (|159[) 
monotonically decreases with c\ as well (cf. Fig. 2]), the eventual vanishing of £ prevents the variational solution 
to attain c\ values too close to unity where our O (ci) treatment of the AU interactions would break down. Hence 
the bounded existence range of £ (fj,, c\) > solutions, inside the critical line of Fig. [5J is indispensable both for our 
approximations to remain valid and for the quantitatively reasonable size prediction of the gluon condensate (|190[) 
and other vacuum scales to emerge. 

Another issue meriting discussion is the order of the transition between the strongly- and weakly-coupled vacuum 
phases. Above we found the (dis)order parameter £ to vanish continuously at the phase boundary, which indicates 
a second-order transition. The same order was encountered in the Ref. [17j ]. although lattice simulations predict a 
first-order transition [43| . In Ref. [l7j this mismatch was argued to be a shortcoming of the mean-field approximation 
(which should break down close to the transition). Nevertheless, one may wonder whether the mean- field approxima- 
tion could generate a first-order transition when the higher- derivative interactions corresponding to c„ ^ are taken 
into account. One could imagine, for example, that the lowest-energy gap-equation solution £ > may cease to exist 
for increasing /x before reaching zero and/or that a competing solution £ ~ may become energetically favorable 
before any of them vanishes. This raises the question whether more than one solution of the gap equation (|119|) could 
exist. In principle, this is indeed possible. An example can be constructed by treating the dominant 2U part of the 
Ci =/= interactions perturbatively to O (ci), instead of resumming it as we did in Eq. (|74[) . In this case two solution 
branches of the gap equation would indeed emerge. One of them is not continuous in the c\ = limit, though, and 
the energy competition with the other one turns out to be unable to generate a first-order transition. In our case such 
scenarios are excluded from the outset, furthermore, since the nontrivial solution of our gap-equation is unique. 



E. Impact of the higher-gradient corrections 

From Sec. [V] onward, we have restricted our practical calculations to contributions from the dominant higher- 
gradient interactions, associated with the low-momentum coupling c\. Hence a few comments on the expected impact 
of the subleading higher-gradient corrections (governed by the c n >2 couplings in the soft-mode Lagrangian (|43|1 ) may 
be useful. All c„ dependence originates from the expansion Q13[) of the infrared covariance. Through the generating 
functional (|45[1 it then enters the soft-mode matrix elements (cf. e.g. Eq. (|103p ) and thereby any amplitude, both in 
terms of the propagators A -1 and A" 1 (cf. Eqs. (|75|) and (| 1 00[1 ) and by virtue of the perturbative expansion l|58p 
of the 4J7 contributions. In order to take the contributions from the 2 (n + l)-gradient interactions with n > 2 into 
account, one therefore has (i) to include the explicit c n dependence from Eq. (|20[) in the expressions for the chromo- 
electric (| 1 37[) and -magnetic (|139[) expectation values, (ii) to continue the expansion of the 4U contributions (j58"| to 
higher c„ , and (iii) to generalize the integrals (|B1[) and (|B2[) , which summarize the impact of the c n contributions on 
the 2U part of the soft-mode amplitudes, to 

!»({,<*,<*,*,...) :=y w 2 (1 _ ClK 2 + C3K 4_ C8K 6 + ... ) + g2 d " n* 1 ' ( 198 ) 



Jn(f,ci,C2,C3,...) := / — — -^dK n > 1. (199) 

J0 [K 2 (1 - ClK 2 + C 2 K 4 - C3K 6 + ...) + £ 2 ] 

(The monotonicity properties of the integrals (|198l) . (|199p are straightforward generalizations of those discussed in 
App.B) 

In Eqs. (|198[) and (|199[) the low-momentum constants c n appear multiplied by 2n powers of the dimensionless 
momentum k £ [0, 1] . This provides the basis for a consistent power counting. Indeed, the contributions from larger 
n are systematically suppressed and at least parametrically dominated by the four-derivative interactions associated 
with c\. Since the c„ are furthermore bounded by the normalizability constraints (|22p and since our optimal value 
c* ~ 0.15 keeps already the dominant contribution rather small, one may expect that a variational treatment of the 
c n >2 contributions would result in even substantially smaller corrections. (Vacuum energy, gluon condensate and other 
matrix elements could therefore be estimated by expanding the integrals (|198p and (|199p around c„> 2 = 0, incidentally, 
with the expansion coefficients determined by the integrals given in App. |Bj) With the above considerations in mind, 
we have not pursued the quantitative evaluation of c„>2 contributions in this paper, although our framework is fully 
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equipped to take them into account. Another reason for neglecting these corrections was that their inclusion would 
reduce the transparency of our explorative study. (The visualization of the vacuum energy density and its minimum 
in Fig. SI for instance, is possible only if c„> 2 = 0.) 



F. Further physical implications and applications 

We conclude our analysis by briefly reviewing several additional features of the resulting vacuum description. To 
begin with the UV physics, we note that the one-loop Yang-Mills j3 function can be recovered almost completely from 
the dynamics implemented in the UV covariance (jlll) (up to a missing 6% correction due to color screening 35] ) . As 
already mentioned, the screening contribution may be accounted for as well, e.g. by generalizing the tensor structure 
of the UV covariance [lH . (We note in passing that Rayleigh-Schrodinger perturbation theory, even to higher orders 
of the gauge coupling, could be set up directly in our wave-functional basis ([TJ) since it renders Gauss' law manifest.) 

Turning to the IR physics content, we start by emphasizing that dimensional transmutation emerges naturally from 
gauge- invariant wave functionals ([1]) with Gaussian cores ©, and that it generates the infrared mass scale and the 
crucial mass gap dynamically. Since a Gaussian becomes the exact ground state of both massless and massive vector 
fields in the non-interacting limit, furthermore, one may expect it to yield an adequate description whenever the main 
net effect of the interactions is to create a quasi-particle-like selfenergy. Similarly, it was argued in Ref. [ItJ that the 
Gaussian approximation should provide a reasonable description of Yang-Mills vacuum physics as long as the latter 
is dominated by a single gluon condensate (which QCD sum-rule analyses suggest especially in the glueball sector 

SHI). 

An important advantage of the gauge-projected wave functionals ([T]) is that they incorporate the nontrivial topology 
of the gauge group and thereby of the gauge fields. In particular^ they contain instanton effects [12, HH which play 
an especially important role in the spin-0 glueball sector [36|, Hg, l50j | and emerge even in basic holographic duals 
for Yang-Mills theory {51] . Furthermore, the functional (J]) comprises large classes of additional, gauge-invariant 
IR degrees of freedom which are often of topological origin as well (2^|. Those are typically collective excitations 
gathered by gauge orbits of vacuum fields which were found to dominate the saddle-point expansion of the generating 
soft-mode functional P5|) and to provide gauge-invariant physical interpretations e.g. for merons and Faddeev-Niemi 
knots [HIH. 

Another crucial consequence of the infrared fields populating the Yang-Mills vacuum, confinement, is characterized 
by the emergence of a linear potential between sufficiently far separated, static quark sources. Despite the rather 
direct access to the interquark potential provided in the Schrodinger representation [IH [20j j and heuristic arguments 
given m Refs. [13, [l8l[20j. however, it remains to be shown that gauge-averaged wave functionals ([I} with a Gaussian 
core of the type are capable of generating the confinement-induced area law for large, spacelike Wilson loops. 
Several cautionary lessons from 3-dimensional, compact electrodynamics can be found in Ref. ;2lj. (The behavior of 
the adjoint Wilson loop would be of interest as well.) Explorations of confinement in Coulomb gauge encouragingly 
found a Gaussian wave functional, divided by the square root of the Faddeev-Popov determinant, to generate a linear 
heavy-quark potential 49] . 

In this context, one should further emphasize that confinement (as well as other crucial vacuum features) is expected 
to emerge in our framework only after gauge projection. In fact, Feynman had argued almost three decades ago Q 
that for a mass gap to be generated in a non-Abelian gauge theory, all field configurations with non- vanishing support 
at "large" A should either be gauge-equivalent to others with support only at "small" A or damped by their magnetic 
field energy. This argument is suggestive because then, in analogy with finite-dimensional quantum mechanics, the 
wave functionals of the ground and first-excited states cannot arbitrarily reduce their energy difference by favoring 
long-wavelength potentials A extending over all of space. As a consequence, the energy or (in the static case) mass 
gap cannot be closed. It would be important to check whether a gauge-projected Gaussian wave functional can realize 
this mechanism. 

In order to simplify the analysis and to focus on the still rather poorly understood infrared gluon sector of QCD, 
we have restricted our investigation to pure Yang-Mills theory without quarks. As the quenched approximation 
to lattice QCD, this provides stable and unmixed glueball states which set a benchmark for sorting out the QCD 
glueball spectrum. (For a first attempt to study glueball states on the basis of a gauge-invariant vacuum wave 



functional see Ref. 521] .1 However, static quark sources could be straightforwardly implemented into the gauge- 
invariant vacuum functional [20l [2l| and even the inclusion of light, dynamical quarks adds no conceptual problems, 
although a nonperturbative treatment of their interactions with the gluon secto r ap pears challenging. (For a variational 
study of fermions in the 1+1 dimensional Abelian Schwinger model see Ref. [531].) 

We close this section by briefly discussing the improvement potential of our approach and by mentioning a few 
additional applications. The 0(g) treatment of the energy density and the O (ci) treatment of the AU interactions 
could be systematically refined. In addition, one may evaluate the corrections due to higher-gradient contributions 
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with couplings c„>2 in the effective soft- mode action (|43|) . For the gauge group SU(2), the mean- field analysis of 
the integral over the auxiliary field £ could alternatively be performed in quaternionic variables, as suggested in Ref. 
[l8j . If this formulation could be generalized to higher N c , it would keep the fluctuations around the mean field under 
control, and the corresponding saddle-point approximation may in fact become exact in the large- N c limit 45] . One 
could further generalize the trial functional basis by including both longitudinal and transverse mode contributions 
(potentially non-analytic) to the covariance, or even allow for a non-diagonal color structure. Finally, one may consider 
to solve the nonlinear soft-mode dynamics (|41[) and its phase structure exactly on a lattice. 

A particularly interesting class of Yang-Mills amplitudes are (equal-time) glueball correlation functions. They 
contain information on the whole spectrum and the decay constants of glueballs which was first explored variationally 
in Ref. [52| and is accessible in our framework as well. It could e.g. be analyzed by comparison with results from 
lattice simulations [54j . the operator product expansion J4(| HH and holographic strong-coupling duals 51]. The 
pseudoscalar glueball correlator contains the topological susceptibility (in the zero- momentum limit), furthermore, 
whose evaluation would complement the analysis of the topological vacuum structure mentioned above. Another 
interesting type of amplitude which may become at least approximately accessible are the Bethe-Salpeter amplitudes 
of glueballs for which e.g. lattice [H3] and instanton liquid model {E3j results are available. A further natural and 
interesting extension of our analysis would be its generalization to finite temperature, as pioneered in Ref. [55| for the 
minimal trial- functional family of Ref. 1171 . This would open up several important applications, including the study 
of the deconfinement phase transition [56| as well as of properties of the currently intensely debated and potentially 
strongly coupled gluon plasma at temperatures of up to 2-3 times the critical temperature [5?} ■ In the longer run, 
the potential of our Minkowski space formulation to describe quantum nonequilibrium physics should be exploited, 
too, e.g. by applying it to currently much discussed nonequilibrium processes in the early universe as well as in the 
aftermath of ultrahigh-energy nuclear collisions (for a brief introduction see Sec. V of Ref. [9J). Hopefully, some of 
the above applications may also provide guidance for how to improve the Gaussian approximation ([3]) to the core 
wave functional. 



IX. SUMMARY AND CONCLUSIONS 

We have implemented and studied a gauge-invariant variational approximation scheme for the Yang-Mills vacuum 
wave functional. Our approach is based on minimizing the Yang-Mills Hamiltonian in a gauge-projected Gaussian 
trial-functional space that allows for an essentially general momentum distribution of the soft vacuum fields. This 
amounts to a substantial improvement upon previously used trial states which populated all infrared modes with 
identical, i.e. momentum-independent strength. Our extension of the trial basis rests on a controlled expansion of 
the infrared- mode dispersion relation (or "covariance"), which characterizes the vacuum wave functional, into powers 
of the soft momenta divided by the dynamically generated mass scale. The momentum dependence can then be 
systematically approximated by truncating the series, which allows us to restrict the quantitative part of our study 
to the leading-order corrections. After complementing the generalized soft vacuum-mode population by the leading 
perturbative ultraviolet modes, finally, one ends up with the largest gauge-invariant trial space for the Yang-Mills 
vacuum wave functional (in 3+1 dimensions) which is currently available and which keeps the variational analysis 
analytically manageable. 

The multifaceted vacuum physics emerging from energy minimization in this extended trial space turns out to 
be partly shaped by a subtle interplay between different parts of the dynamics. To begin with, the leading-order 
momentum distribution of the infrared modes is characterized by just one variational parameter. After translating 
the dynamical content of the wave functional into an effective Lagrangian for collective sets of soft gauge-field orbits, 
this parameter reappears as the coupling constant of the dominant higher-gradient interactions which modify the 
low-momentum dispersion of the vacuum modes. Moreover, it acquires an intuitive physical interpretation in terms of 
the group velocity of gauge-field quasiparticles in the vacuum. While this velocity vanishes in the unimproved ground 
state, one of our main findings is that it becomes finite in the improved one. The novel dispersion relation further 
reveals that the infrared quasigluons in the vacuum decelerate when a force is applied to them. In other words, their 
"effective kinetic mass" is negative (analogous situations are encountered in several condensed matter systems) and 
their color conductivity can increase with increasing strength of an external color-electric field, implying a negative 
differential resistance of the vacuum. 

The variationally optimized value of the higher-gradient coupling turns out to be positive. This sign may at first 
appear surprising since it is opposite to the one induced by the massive vector propagator which was previously 
adopted as a model for the infrared covariance. Its physical impact is quite intuitive, however, since it implies 
that soft gauge-field modes with larger momenta populate the vacuum more strongly than their longest-wavelength 
counterparts. Vacuum fields with constant color-magnetic field strength are thus energetically disfavored, in particular, 
which prevents the onset of the Savvidy instability. An energy barrier against larger quasigluon velocities preserves 
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the normalizability of the vacuum wave functional and keeps the variationally determined value of the quasigluon 
velocity small. In addition, it implies the existence of a saddle-point solution for the disorder-parameter field at a 
dynamically generated mass scale which is consistent with lattice results for the Yang-Mills scale. (As a practical side 
benefit, the relatively small optimized coupling-parameter value also justifies a leading-order pcrturbative treatment 
of the four-soft-mode interactions.) 

Moreover, the above findings indicate that the dimensional transmutation mechanism is an almost quantitatively 
robust feature of gauge-projected Gaussian core functionals (instead of emerging accidentally from an oversimplified 
trial space). As a consequence of this rather nontrivial result the lowest-dimensional gluon condensate, which scales 
with the fourth power of the dynamically generated infrared mass scale, acquires a magnitude well inside the range 
predicted by other sources. This observation further supports the conjecture that Gaussian trial functionals provide 
a reasonable vacuum description in dynamical situations where just one vacuum condensate dominates. While this 
holds e.g. for the Cooper-pair condensate in the BCS theory of superconductivity, it is less well established in Yang- 
Mills theory where the lowest-dimensional gluon condensate provides the natural candidate. The main reasons for the 
robustness of the vacuum scales turn out to be that the nontrivial solution of the gap equation is unique and that the 
disordered vacuum phase exists only in a limited domain of all variational parameters. Even the soft-mode dynamics' 
order-disorder phase transition, where the vacuum energy attains its minimum, remains squarely of second order. 
This is in contrast to the first order predicted by lattice simulations and probably a shortcoming of the employed 
mean-field approximation near the transition. In any case, the improved dispersion generates a considerably larger 
attraction among the variationally optimized vacuum fields and lowers the vacuum energy density by more than 10%. 
This indicates that the nonvanishing quasigluon velocity provides a rather substantial overall improvement of the 
vacuum description. 

Nevertheless, the gauge-invariant variational framework remains analytically tractable, the computational costs 
are comparable to those of fully gauge-fixed approaches, and the rather high degree of transparency owed to the 
explicit representation of the vacuum state is maintained. In combination, these qualities suggest considerable further 
potential of our approach as an intuitive theoretical laboratory for gaining insights into the nonperturbative real-time 
physics of Yang-Mills theories. The established framework provides, in particular, a privileged testing ground for the 
impact of topological vacuum excitations and their relation to the gauge symmetry. Indeed, both aspects are manifest 
in our effective soft-mode dynamics and should be explored further. In addition, our approach seems well suited to 
chart the maximal physics content which gauge-projected Gaussian core functionals can accommodate. In the longer 
run, this would allow to clarify the principal limitations of the Gaussian approximation and may help to trigger new 
ideas for going beyond it in a both systematic and analytically manageable way. 
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Appendix A: Evaluation of vacuum expectation values 



In the following appendix we review the general approach to calculating vacuum matrix elements in gauge-projected 
trial functionals of the type (HJ, following Ref. [l7|. We recast this framework into the language of generating 
functionals and also derive several expressions which will be needed in the main text. 

To begin with, consider local gauge-invariant operators O (A, E) composed of static gauge fields Af (x) and their 
canonically conjugate momenta —Ef (x) (where E is the chromoelectric field) in the Hamiltonian formulation of 
SU(-ZV C ) Yang-Mills theory, adopting temporal (or Weyl) gauge. In the Schrodinger coordinate representation, the 
canonical commutation relation [Ef (x) ,A^ (y)] = iS ah 5ijS 3 (x — y) then implies that Ef (x) — iS/5Af (x). Our aim 
is to calculate the (trial) vacuum expectation value of O, i.e. the matrix element 



(0(A,E)) = 



JDA% 


A 






A 


JDA% 


A 




A 





(Al) 



where the unrestricted linear measure DA comprises all gauge orbits and topological charges. (Recall that Gauss' law 
is enforced as a constraint on the physical states and in particular on the vacuum wave functional vPq, cf. Sec. Ill AH 
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1. Exact integration over the gauge field 

After inserting the (residual) gauge- invariant, unnormalized wave functional QJ into Eq. (jAlll (and interchanging 
the order of the A, U and U integrations), the matrix element becomes 



(0(A,E)) = 



JDU JDUJ DAipo 


AP 


O (I 


a it 
' SA 






jDUjDUjDAtpo 











(A2) 



Since O is gauge-invariant, the integrand in the numerator can only depend on the relative transformation U := U U. 
Indeed, after changing the integration variable A — > A' — A u (for fixed U), one has 



DAip 



A' 



A' 



= / DA'i; 



A 



iu- 1 u 



0(A',^)^A'\ 



(A3) 



and the analogous expression for the denominator. Making use of O (A) = O (A') and DA = DA' , renaming A' to A 
and employing the translational invariance DU — DU of the Haar measure then yields 



DU / DAiPo 



A 



A 



= DU DAi(j 



A" 



(A4) 



(and the analogous expression for the denominator) which does not depend on U since the latter was absorbed into 
the integration variable A. Hence the factored (infinite) gauge group volumes J DU in numerator and denominator 
cancel, leaving us with 



JDUJ DAiP 


'a u 




■>Po 


A 


JDU J DAiP 


A u 


■00 


A 





(0(A,E)) = 



Now it is useful to define an effective bare action r& [U] and the associated nonlocal gauge-field correlators as 



(A5) 



exp{-r b [£/]} := / DAiP 
A...A...E...E)))ex V {-T h [[/]}:= / dA% 



A 



A 



=: I DAe-^ u \ 
A , 



7 ? id iS , 
A...A...— ip 
SA 6A 



(A6) 
(A7) 



where the triple-bracket (((...))) denotes (nondiagonal) matrix elements between [/-rotated and unrotated core- 
functional states. (Note the symmetry 1\ [U] = I\ g^Ugn with gh.n G SU(iV c ) which originates from the transla- 
tional invariance of the two group integrations in Eq. (| A2[l ■) With the relation 



O A\ 



iS 



6A a (x) 



tpo 



A 



O [A] V-o 



.4 



(A8) 



which replaces the (quasi-) local, differential operator O by a nonlocal functional O of the A fields only, one then 
arrives at 



(0(A,E)) 



JDU J DAipo 




O [A] V> 


A 


fDU{{(0(A,E))))exp{-r b [U]} 


JDU J DAipo 


a u ~ 




A 


J DUexp{-T h [U]} 



(A9) 



Hence the calculation of the matrix element is reduced to evaluating the vacuum expectation value of the expressions 
(|A7[) . specialized to contain A fields only, in the U field dynamics determined by 1\. 

Owing to the Gaussian nature of the wave functional (|3]), the effective bare action defined in Eq. (|A6[) and hence 
the matrix elements (|A7[) can be calculated exactly. In order to render the U dependence explicit, we write 



j[A,U] :=-]*% 



A 



4>a 



A 



1 r 



(A u y G- 1 A U + A T G- 1 A 



(A10) 
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(in a triple matrix notation for the color, vector and spacial "indices") where the rotated gauge field 



(l* 7 ) a = S ab (x) A b + g- x L a (x) 



is expressed in terms of the homogeneous and inhomogeneous transformation functions 

S ab (x) = ifr {A a [/ t (x) \ b U (:r)} , L a (x) = itr {a !/ 1 (x) dU (£)} . 



(All) 



(A12) 



After inserting the representation (|A11|) into Eq. ()A10j) . one finds that 7 contains an A independent piece and a term 
which is bilinear in A, 



Above we have defined 11 



7 [A, U] = - (A + af M (A + a) + -^L T VL. 



M := S T G- 1 S + G -1 = M T , a := g- 1 M~ 1 S T G- 1 L, 
V := G- 1 - G- 1 SM' 1 S T G- 1 = (G + SGS' 1 )' 1 . 



(A13) 



(A14) 
(A15) 



(G 1 and M. are symmetric in a, 6, i,j and x, y, and transposition of a is understood to invert the direction of the 
derivatives in G _1 as well.) After performing the Gaussian integration over A, we arrive at the explicit form 



r b [U] = / d 3 x / d 3 y£ h (x, y) 



In / DAvq>{- 1 [A,U]} = -Tr]n(M) + —s 



L T VL 



(A16) 



of the bare action which is a nonlocal, nonlinear cr-model. (The regularization by the ultraviolet cutoff Ayv is implicit 
in the definition of the wave functional.) This dynamics sums up the self-interaction of the "background" field U 
which the gauge fields induce through the vector a[U] and the propagator M^ 1 [U]. In the following we will omit the 
first term oc Tr In (A4) in Eq. (|A16|) since it is of O (<? 2 ) relative to the second one and does not contain U derivatives 

The connected (equal-time) n-point functions at fixed U, 



(((A%(z 1 )A%(z 2 )...A? : (z n )))) = 



J DAiPv 


AF 






..A 




A 


JDM 




■00 


A 





(A17) 



specialize the definition (|A7[) exclusively to A fields and can now be evaluated by means of the generating functional 



DAiPo 



A 



A 



exp(-JA) = / Z?Aexp(-7 



1 



A,U 



- J A 



exp ( -JM~ 1 J + Ja ) exp (-T h [U]) 



as 



{{{A%{z x )A%{z 2 )...A%{z n )))) = 



(-1)" 5 n exp U JM- 1 J + Ja 



5J£ (fi) ...8Jt: (z n ) 



(A18) 
(A19) 

(A20) 



j=o 



This expression reproduces the Wick expansion for correlators of a vector field A with nontrivial propagator A4 _1 
in the background of the U and a fields. (In contrast to free or mean-field correlators, it does not factorize 2n-point 
functions into products of n 2-point functions, however, and also generates finite (2n + l)-point functions.) Explicit 
expressions for the 2, 3 and 4-point functions needed in the main text are 



(((ah (*o a z m)) = x^r 2 20 + <i &) <i &) 



(A21) 
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and (in shorthand notation for arguments and indices) 

(((A (1) A (2) A (3)))) = - [M- 1 (2, 3) a (1) + Mr 1 (1,3) a (2) + M' 1 (1, 2) a (3) + a (1) a (2) a (3)] (A22) 
as well as 

(((A (1) A (2) A (3) A (4)))) =X- 1 (1,2)X- 1 (3,4)+M- 1 (1,3)X- 1 (2,4)+M- 1 (1,4)X- 1 (2,3) 

+ AT 1 (1, 2) a (3) a (4) + AT 1 (2, 3) a (1) a (4) + M~ x (3, 4) a (1) a (2) 
+ .A/T 1 (1,4) a (2) a (3) + .A/T 1 (1, 3) a (2) a (4) + X" 1 (2, 4) a (1) a (3) 
+ a (1) a (2) a (3) a (4). (A23) 

Hence the above results reduce the calculation of the matrix element (IA2[) to integrating over the nonlocally interacting 
U (x) fields according to Eq. (|A9j) . This integral will be performed in the next two sections. 



2. Perturbative integration over the hard U modes 

Since we are interested in amplitudes with soft external momenta \pi\ <C fi (where fi is a typical hadronic scale, 
e.g. the lowest glueball mass), the integrations over U in Eq. (|A9|) can be approximately done in two steps [17|: after 
integrating out the hard modes U> (containing momenta k > /i) perturbatively, one performs the integral over the 
remaining soft modes U < (which contain the momenta k < fi) in the saddle-point approximation. In preparation for 
this procedure, we factor the integration variable U in a group-structure preserving manner as 

U (x) = [/< (x) U> (x) , U> (x) = exp {-ig<p a (x) A Q /2) . (A24) 

Note that the soft modes vary only weakly over distances \x — y\ < fi , i.e. U < (x) ~ U < (y). Since asymptotic 
freedom is manifest in the high- momentum covariancc (jllj) . the bare coupling gj, is small at the large cutoff scale 
Auv 3> Aym where the theory is defined. Hence the hard modes can be integrated perturbatively e.g. by Wilson's 
momentum-shell technique [581 ] , down to the infrared scale /i chosen such that the renormalized coupling g (/x) remains 
sufficiently small. The functional measure factorizes according to Eq. (IA24[) as 

DU = DU < DU > oc DU < D4) a . (A25) 

As a consequence, the bare action becomes a functional of the <j) and U < fields, 

r b [</», [/<] - 1\> [0] + r b ,<> [</>, + r b .< [u<] , (A26) 

and Eq. (|A9|) for the matrix element turns into 

tn(A I DU < SJ^WP (A E)))) exp {-r b [0, U K }} 

As long as <7 2 (/«) <C 1, the field can be integrated perturbatively to one loop in the bare coupling (where the integrals 
over loop momenta are regularized by the cutoff Auv)- The result has the form 



((O (A, E))) exp {-r< [£/<]} := / (((O (A, £)») exp {-r b [0, {/<]} (A28) 

where we defined the soft-mode action T< [t/<] as 

exp {-r< [[/<]}:= y^cxp{-r b [0, [/<]}, (A29) 

made use of r bj<> [J7<, 0] = to leading order in the coupling, and assumed that the (one-loop) anomalous dimension 
of O vanishes. (This is appropriate for the conserved Yang-Mills Hamiltonian density (|131[) which we consider in the 
bulk of the paper.) The expression (|A27j) for the matrix element is then further reduced to the soft-mode integral 

iO (A F)\ JDU<((Q(AE)))C XP {-T < [[/<]} 

{ ° {A E)) ~ jDU < e^{-r < [[/<]} ■ (A30) 
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In order to exhibit the (approximate) U < dependence of T< [t/<] and ((O (A,E))), we now exploit the smallness of 
the bare coupling to truncate the expansion of /7> in powers of g. Up to corrections of O (g 3 ) one finds 

u> = ex P (-igr^r) = i - - C W>° Aa ) 2 + o (g 3 ) , (A31) 



a 2 

t = di<j) a + O (g 3 ) , S^ = 5 ab +gf abc <p c -i-f ace f bde cl) c cl) d + 0(g 3 ), (A32) 

A? - S« 6 A* , 4 + X a >ti = X a <A + d^ a + gf^X^ - ^■f ace f bde X h < i ((> c (t> d + O {g 3 ) . (A33) 

(Above we have rescaled the Cartan-Maurer form (|35|) as Xf := Lf/g in order to facilitate explicit power counting 
in g. The spacial index should avoid confusion with the SU(iV c ) generators.) To the same order, the operator (IA15|) 
expands as 

vab ( ~ g) = G-i(x-ff) sab 1 1 + _g^_ [r (f) ^ (f) _ ^ (f } ^ a ^ + ^ ^ m | 

+ O (g 3 ) . (A34) 

(Note that V is independent of [/<, as a consequence of S^GS^ 1 ~ G which holds because 5< varies slowly over 
distances |x — y| < over which the decaying G remains non- negligible.) Inserting the above expansions into the 
bare action (IA16|1 . one finds the associated Lagrangian 

Lb (x, y) = Cb,> (x, y) + £ b ,< (x, y) + £ b ,<> (x, y) + O (g 3 ) (A35) 

with the bilinear parts 

£b,> (S, y) = \di<P a ($) G" 1 (x-y) (y) , (A36) 

£ b ,< (x, y) = \\% ti (x) G- 1 (x-y) X a <A (y) (A37) 

and the hard-soft-mode interactions (after using G _1 (x — y) = G _1 (y — x) and anticipating the exchange symmetry 
x «-» y after integration over x and y) 

A,<> {S, y) = 9 -f abc X b <ti (x) G- 1 (x - y) <f> c (x) d^ a (y) 

g 2 { N c 6 ab 6<* bde 
+ 8 ^ iV 2 - 1 11 
x X b ^ (x) G- 1 (x - y) X a <ti (y) [<f> c (x) 4> d (x) - cf> c (x) <j> d (y)j . (A38) 

(Note that £ b is invariant under <p —¥ —<fr.) To leading order in g the hard and soft modes decouple, and the higher- 
order interactions £ b ,<> can be treated perturbatively. (The O (g) term in £ b ,<> contributes neither to the /3-function 
nor to the soft-mode action.) Writing further 

r b [u] = r b ,< p7<] + r b ,> [</>] + g v 1 [<f>, u K ] + , 9 2 r 2 [</>, u<] + o ( g 3 ) (A39) 

renders the bare-coupling dependence of the action explicit. In terms of the kinetic operator 

K ab (x - y) = ±6 ab [dyAtG- 1 (x - y)} p2>/j2 (A40) 
of the Lagrangian (|A36j) . furthermore, the hard- mode action becomes 

T b ,> [0] = J d 3 xj d 3 yC >>b (x , y) = ±jd 3 xj d 3 yr (x) K ab (x ~ y) ^ (y) + O (g) (A41) 

= 3 (^ [e (A&V " fc2) " 6 fc2)] ~ r ( fc T (fc) k2 ~ r (~ % ) (A42) 



3G 



from which one reads off the static (j) a propagator 



J r>0exp{-r bj > [4>]} 

d3k [6 (A 2 VV ~k 2 )~9 ( M 2 - k 2 )} j^-v) 2G >j k h ab 



(27T) 



k 2 



(A43) 
(A44) 



The nonstandard kinetic term of the high-momentum a model (|A42[) thus generates a propagator with large-A: asymp- 
totics K^ 1 (k) oc G> (k) jk 2 ~ 1/fc 3 (instead of the standard 1/fc 2 behavior which impedes perturbative renormal- 
izability in three dimensions). Hence the high-momentum dynamics (IA42I) is renormalizable (the tadpole diagram 
diverges only logarithmically) and asymptotically free, with a perturbative ft function very similar to that of Yang-Mills 
theory 



From the definition (|A29I) and the decomposition (|A35[) of the bare action one then has 

exp{-r< [[/<]} = exp{-r b ,< [[/<]} j £>0exp{-r b ,> [</>] ~ r b ,<> (A45) 
= exp{-r b>< [[/<] - (r b ,<> [C/<,0]) rb> +0(g 2 )) (A46) 
(where the average (••■)r b > over <f> is weighted by exp (— r b > [<fi])) which yields the renormalized soft-mode action 



r< [u<] = -r b>< [c/<] - (r b ,<> [0, t/<] 



O (g 2 



l 



Ag 2 ( M ) 



d 3 x / d 3 yL a <ti (x) G- 1 (f - y) (y) + O (g 2 ) 



(A47) 
(A48) 



To the considered order (and for G> (fc) = k 1 ), T< has therefore the same form as the bare soft-mode action I\<, 
but with the bare coupling replaced by the renormalized one, 



9 O) = 9b + 



glN c j Auv 



(27T) 



(A49) 



which was calculated in Ref. 35| . This renormalization of the soft-mode interactions compensates for the removal of 
the hard E/> modes from amplitudes with external momenta p 2 -C /J 2 ■ Analogously, the <\> mode contributions 



((0(A,E)}) 



J Dcj) ({{O (A, E)))) exp {-r b [[/<, <£]} 



jD^>exp{-T b [[/<,</)]} 



can be integrated out perturbatively by rewriting 

exp {-r b ,< [£/<]} 



({0(A,E)}} 



exp{-r< [[/<]} 
jD<t>(((Q (A, E)))) exp {-r b ,> [</>]} 
/D0exp{-r b!> [0]} 



(A^))))exp{-r b , <> [[/<,</>]} exp {-r b ,> [</>]} 



(A50) 

(A51) 
(A52) 



(To the considered order the anomalous dimension of O {A, E) arising from the perturbative interactions T b <> [U < , cj>] 
can be neglected. This approximation becomes exact for the operator considered in the main text, i.e. the conserved 
Yang-Mills Hamiltonian density (|131l) in temporal gauge.) After introducing another generating functional 



z> [j] ■■= J £>0exp 
= z y [0] exp 



-r b ,> [cj>] - J d 3 xj a (x) r (x) 



(A53) 
(A54) 



the <fi a fields in (((0 (A,E)))) can be replaced by derivatives —5/5j a with respect to the source j a , which yields 

1 



((0(A,E))) 



[0] 



z> [j] 



(A55) 



and shows that this intermediate matrix element is a functional of G 1 and U < only. 
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3. Saddle-point integration over the soft U modes 

With the result (|A55[) for the intermediary matrix element ((O (A, E))) at hand, it remains to perform the functional 
integral 

J £'t/<cxp{-r < [Z7<]} 

over the soft modes U < . In order to prepare for this last step, we integrate over a complex matrix field V and represent 
the unitarity constraint J7<£/< = 1 by inserting a delta functional. The latter is rendered explicit by the additional 
integration over an auxiliary Hermitian field E (x) which acts as a Lagrange multiplier, as explained in Sec. IIV Al 
The result is 

(o (A, E)) = I^{ DV f (^)))ex P { - r< [y] 

where the integration contour for £ (x) runs parallel to the imaginary axis and where 

T E [V, n = ^J d 3 xtr {E (V*V - 1) } . (A58) 

(Note that E is Hermitian since V is. The unimodularity constraint det £/< = 1 could additionally be implemented 
by integrating over a minimally coupled U(l) gauge field [l^. We refrain from doing so because the impact of the 
difference between SU(A r c ) and U(N C ) is of order 1/N 2 .) The matrix element (|A57I) can then be obtained from the 
soft-mode generating functional (jT7| . i.e. 



Z[j,?] = J DzJ DVcxp -r< [V] - T s [V, E] - / rf 3 ztr {jyt +jty} 
after replacing the J7< and fields in ((O (A, S)}} by derivatives with respect to the matrix sources and 



(A59) 



(A60) 



The unconstrained integral over V in Eq. (|A59[) is Gaussian and can be done exactly (cf. Sec. IIV Al) . and the 
remaining integral over E can be performed in the saddle-point (or mean-field) approximation, as detailed in Sec. 
IIV Bl To leading order, this just amounts to substituting the solution E of the gap equation (|62|) into the integrand 
ofEq. (|A59|> . 

Given an explicit expression for the covariance, as in our case Eq. ((T3"|) . the soft-mode action T< [£/<] (cf. Eq. (|4Tj) ) 
and the amplitudes (|A60j) are thus uniquely determined functionals of G^ 1 . In the (trial) vacuum expectation value 
(Hym) oi the Yang-Mills Hamiltonian, in particular, G^ 1 plays the role of a variational trial function. By means of 
the parametrization (|20j) . ("Hym) becomes a function of the variational parameters /x, m g , and {ci} with respect to 
which it has to be minimized according to the Rayleigh-Ritz procedure. (The bare cutoff Auv 3> Aym, on the other 
hand, has been traded for Aym which will be fixed from lattice data, cf. Sec. I Villi ) 



Appendix B: Coincidence limit integrals 

The coincidence limit of the soft-mode correlators calculated in Sees. IIVI and Ivl is (for c n >2 — 0) determined by 
integrals of the type 

U^)^ KHl _ ClK 2 )+e ^ n>l, (Bl) 

3n — dn n>l (B2) 

JO [K Z (1 - ClK Z ) + ^\ 



with c\ < 1 (which ensures normalizability of the vacuum wave functional). In the following we list pertinent properties 
of these integrals and derive analytic expressions to be used e.g. for the numerical solution of the gap equation and 
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for the evaluation of the vacuum energy. To start with, we note the obvious relations 

H (£.ci) > i n +i (£,ci) , (B3) 
in (f> ci) > j n+ i (£, ci) (B4) 

which hold point by point. Furthermore, for the physically required ci < 1 all integrals (|Bip . (|B2|) are non-negative. 
They decrease monotonically with increasing £ > and increase monotonically with increasing c\ in the allowed range 
— oo < ci < 1 (where 1 — c\k 2 > decreases with increasing c\ for any k G [0, 1]). All integrals are regular at c\ = 
and at £ = 0, with the exception of ji which contains a single pole at £ = (see below). Due to the factors of n 2 in 
the numerator, all integrands have most of their support close to k — » 1, i.e. at momenta k ~ (i. 



1. Analytical expressions 



An efficient strategy for evaluating the integrals i n (£, ci) analytically starts from the identity 

1 1/11 



K 2 (1 — ClK 2 ) + £ 2 y/l 4- 4ci^ 2 V^ 2 — K 2 K 2 — K 2 

which exhibits poles at the positions 



Hence the evaluation of 

a/1 + 4ciC 2 ^0 V « 2 - K 2 K 2 - K 

is reduced to the calculation of one-pole integrals. The simple identity 



(B5) 



«i, 2 feci) = — (l±Vr+4ci^j. (B6) 



„2n *2(n-l) 

^-- 2 ^ + « 2(n - 1) (B8) 



1 K 2« 2 /-l K 2(n-1) 

T^— = TT-^ + x 2 / -r— o dn (B9) 



implies the recursion relation 

2n- 1 

and allows to reduce the l n (£, ci) to analytic expressions involving only 

Q 1 Q rf« = — -In ^4 = --arctanh fi^ (for a; 2 g 10, (BIO) 

— x 2x x — 1 a; V^/ 

(Note that our n\ 2 ensure a; 2 £ ]0, 1[ inside the k integration range (except at £ = 0).) From the recursion relation 
one then finds (for x 2 g" ]0, 1[) 

1 K 2 / 1 \ 

dn — 1 — a; arctanh ( — I , (Bll) 
\xj 

dn = - + x 2 - x 3 arctanh [ - J , (B12) 



K 2 


— X 2 










K 2 


— X A 




K 6 



3 V x 

1 1 

5 + 3 



drc = ~ + -a; 2 + x 4 - x 5 arctanh (-),•■• (B13) 
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so on. Combining the above results one obtains for the first four i n , which are needed in the main text, the 



and 

following expressions 

*i (C,ci) 
h (£,ci) 
h (?,ci) 



k,\ arctanh {^~^ ~ K 2 arctanh 
v/l+4ci€ 2 
1 — K2 arctanh (-^J 



K 



hi 9 



| + k 2 . — ^2 arctanh (j^j 



1 — «i arctanh (7^") 



| + k 2 — «;f arctanh ( 



(B14) 
(B15) 
(B16) 

^=^= . (B17) 

(Note that the main £ and ci dependence of these expressions originates from the pole positions Ki t 2 (£, ci), cf. Eq. 

CSS) .) 

Analytical solutions for the integrals j„ (£, ci) can either be derived from those for the I n (£, C\) by using the identity 

(B18) 



1 i k: 
5+3 



+ k| — «2 arctanh (;rj) 



| + ^ + «f — k\ arctanh (j!r^j 



1 d 

jn (£,Cl) = ~^-g^ l n (£> C l) ! 



or again by direct integration of the pole decomposition. Following the latter path, we have from the square of relation 
(|B5l) 



(£jCl) 



1 + 4ci^ 2 



which we evaluate further with 



K 2n 2k 2 " K 2n 

> 2 -* 2 ) 2 ~ (« 2 - «?) (« 2 - «i) + ( K 2 - « 2 f 



/■I 2 



1 1 



o (k 2 -x 2 ) 2 2s 2 -1 2x 



-1 arctanh I, 



a; 



/o («a-x?)(«2_aJ) 



a;? arctanh — — ii arctanh — 

_~ £2 £l 



2 2 
X l ~ X 1 



and 



1 3x 2 - 2 3 

rCEK = — ■ 



tui, ; x arctanh — 

/o (k 2 - x 2 ) 2 2 x 2 ~l 2 



/o (k 2 - .x 2 ) ( K 2 - a; 2 ) 



a;? arctanh — — in arctanh — 
d« = 1 - ^ 5 51 



2 2 
•^1 — • Z '2 



as well as 



/o (k 2 — x 2 ) 2 

.,6 



dn = 



2 + 10s 2 - 15x 4 5 , 1 



Do 1 

—a: arctanh — , 

2 x' 



6(1 -a; 2 ) 

1 a;f arctanh ^ — x\ arctanh 

/ /„2 Z T 2\ ( K 2 _ T 2\ dK = o + ( X l + 2 _ T 2 1 

' V 1 / v 2/ 12 



(B19) 

(B20) 
(B21) 

(B22) 
(B23) 

(B24) 
(B25) 



(all for a; 2 ^ ]0, 1[), and so on. 
This yields 



1 



1 + 4ci£ 2 
2k 2 



11 ,11 

7T 

v 2 



2 X 2 1 



2«! 



2 2 



/l^2 2/Cr J 



] arctanh • 



1 



Zft2 1 \ 1 

—5 5- H arctanh — 

K\ — «2 2^2 / K2 



(B26) 
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and 



1 + 4 Cl e 2 
3 

7:^1 - - 



1 3k 2 



2 AC? — 1 



1 3k| - 2 



2k? 



arctanh • 



: K 2 



arctanh ■ 



«2 



(B27) 



and 



h (C,ci) 



1 



1 + 4ci£ 2 



1 15/cf - 10k 2 - 2 1 15k| - 1Q«| - 2 2 k 2 + 3k? - k 2 , - 3k? 
6 k? — 1 6 / 



2k? 



arctanh ■ 



1 



Ki 



- 1 

2/«2 



Ki Kn 



arctanh • 



1 



H 2 



(B28) 



and 



and finally 

Is i€,ci] 



fa (£,ci) 



1 



1 + 4ci£ 2 



16 + 14kJ + 70k? - 105k? 1 6 + 14k| + 70k 



30 



1 — K 2 



15 



2 1 



2k? 



arctanh • 



Ki 



30 



2k^ 



1-K 2 



2 ^2 



arctanh - 



K 2 



105k§ 



(B29) 



1 + 4 Cl £ 2 



1 10 + 18k? + 42k? + 210k? -315k? 1 10 + 18k 2 . + 42k| + 210k| - 315k| 



70 



1 - k{ 



70 



1-K 2 



2 15k? + 21k? + 35k? + 105k? - 15«| - 21k| - 35^ - 105k| 



105 



2 Kl 



2k? 



arctanh 

Kl 



2k 2 



arctanh ■ 



k 2 



(B30) 



(the Ki ; 2 (C> c i) are defined in Eq. (|B6[1 ) which we needed to evaluate the matrix elements derived in the main text. 



2. Limits 



When analyzing the soft-mode contributions to gap equation and energy density, it is useful to consider the £ — > 
and c\ — > limits. This boils down to finding the limits of the i n and j n integrals which we provide in the present 
section. The £ —> limits are (for ci < 1, as before) 



d-K) 1 



? n (f = 0,ci)=/ -dn=- 2 Fi[l,n- -,n+-,ci) (B31) 

jq 1 — C\K ZTl — 1 



In - 1 

for n > 1/2 and 



(B32) 



rX ^.2n — 4 j / 3 1 



ci^O 1 



2n 



(B34) 



for n > 3/2. (The hypergeometric functions 2 -Pi (a, &, c, 2) are defined e.g. in Ref. |59j|.) Note that only ji, which 
does not appear in the coincidence limit of the soft-mode correlation functions, contains a £ —¥ divergence (cf. Eq. 
(|B38|) ). The £ -> limit of the l n integrals can be reexpressed via continued partial integration as 



i \ cy arctanh -,/ci \ /T ™^ 

-(f^)— 5=1 lijjETT ^rH (B35) 
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which yields 



arctanh s [cy c _vo 

h ? = 0, a) = -=J-^ ^ 1, (B36) 



1 / arctanh */c7\ 
F 2 (e = 0, Cl ) = 1 — ^ , (B37) 



ci V \/ci 



etc.. Furthermore, 



as well as 



h {M ^° ~ - ^ - 3,/cI arctanh ,/cl) - ^ Cl ? + O (? 2 ) (B38) 



- . ,1/1 arctanh ,/cT\ 

M£ = 0,ci) = - i + , B39 

2 \1 - ci ^/cl / 

- 1/1 arctanh ,/c7\ _ s 
J3 (? = 0, ci) = — =^ , (B40) 



2ci V 1 - ci v/cl 

1 / 3 — 2ci arctanh v /c7 
2^f ^ 1 - ci ~ 3 V^T 



J4 (? = 0, Cl ) = A ( ^-T 1 - 3 ^ V 1 ) , (B41) 



- , 1 / 15 — 10ci — 2c? arctanh ,/cT\ 



e (?, ci) := l 2 - 2cii 3 + c 2 ?4 + 2ci-z 2 (j 3 - 2cij 4 + Cifa) (B43) 
5=o 1 _ ci _ 27 / _ arctanh ^ \ 



This implies, in particular, 



3 5 3C V yfci 

for the combination of integrals which appears in the electric soft-mode contribution (|149p to the vacuum energy 
density (fl59"]> . 

We further examine the ci — > limits 

f 1 n 2n 

in (?) := In (?, ci = 0) = / d/t, » > 1, (B45) 

ft IS 



1 

o (« 2 +? 2 ) 



Jn (0 := J„ (?, ci = 0) = / ; rf«, n > 1 (B46) 

_1_ 1 



which we distinguish from the general integrals (|B1[) . (IB2|) only by their arguments. Both i n (?) and j n (?) decay 
monotonically in ? £ [0, oo] , starting from the (finite) values 

in(0) = ^—r, n>l, (B47) 



2n- 1' 

i™(°) = ^r-^' ™> 2 ' (B48) 
2n — 3 



except for 



h (?) ^° ~ - 1 + \e - \e + o (? 6 ) (B49) 



which exhibits the already mentioned pole divergence at ? — > 0. The i n (?), j n (?) further satisfy the recursion relations 



-/ di n (?) 



4 (?) = = -2?j n (?) < 0, (B50) 

i;(?) = ^^ = -4?fe n (?)<0, (B51) 
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and 



~ ln (0 = 2^=T ~ (0 ' (B52) 



Jn (0 ' 



(2n-l)e 2 j„-i 



2n - 3 

which may be summed up into a finite power series and a transcendental piece. 

n-2 / tf 2fc 



(n > 1) (B53) 



mo - E Sri + <-"" +1 (' - f 2 - 3 - 



as well as 



(B54) 



k=0 

Explicit expressions for the first few i n (£) and j n (£), which appear in the main text, are thus 

* i (0 = 1 — £ arctan ^ , (B55) 

? 2 (0 = i-C 2 +e 3 arctani (B56) 

h {0 = \~ \ e + ^ ~ e arCtan \ ' (B57) 

* 4 (o = \ - \e + \e -e+e arctan ± (bss) 



ii (0 = ^ (- YT? 2 + 1 arctan ?) ' (B59) 

J2 (0 = i (3 - - 3£ arctan ^ , (B60) 

.73 (0 = \ (| - 5C 2 - + 5£ 3 arctan ^ , (B61) 

.74 (o = ^ Q - \e + Y~re ie arctan ?) ' (B62) 

1/9 3 _ - „ ~ .0 1 



J 5 (0 = 2 y y - 5 (3 - 5£ 2 + 15£ 4 ) e 2 - y~^2 + ^ arctan J ) ■ (B63) 

Finally, we note that the (£), j„ (£) (and their analogs with additional powers of the integrand's denominator) 
determine the coefficients of the expansion of the integrals i n (£; ci) and j n (£; Ci) in powers of ci. For example, from 

1 1 K 4 K 8 

^tt; 57 To = ~5 To + 2 c l + s c l + ( c l) . ( B64 ) 

K 2 (i-ci K 2 )+e 2 « 2 + e 2 (k 2 +£ 2 ) 2 (K 2 +^) 

1 2 "'' + — ^-_c 2 +0( C ?) (B65) 



( K 2 (i- ClK 2 ) + e 2 ) 2 (n 2 + ef (K 2 + ef (« 2 +e 2 ) 

one has (with k n (£) := k 2 "/ (k 2 + £ 2 ) 3 dn) 

In (£ Cl) = In (0 + jn+2 (0 Cl + fc n+4 (0 4 + (c?) , (B66) 

jn K; ci) = j„ (0 + 2fc„+2 (0 ci + O (cl) . (B67) 
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FIG. 1: The generalized IR covariance G< x (fe) = /i (l - cifc 2 //i 2 ) / (1 - ci) for ci = 0.17 (full line) and ci = -0.17 (dotted), 



together with the UV covariance G> x (k) = fc, and in comparison with G 1 (k) — k (medium dashed), G ln g= ^ (k) = \/k 2 + fx 2 
(long dashed) and G^k (&) = M (short dashed) (all curves in units of Aym and for jj, = 8.8Aym). 
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FIG. 2: Vacuum phase diagram: inside the plotted phase boundary p c (ci) /Aym the theory is in its strongly-coupled disordered 
phase. (Our treatment is approximately valid in the range fj, > 4Aym and ci G {—0.5,0.5}). 




FIG. 4: The energy density e(/x, ci) of the vacuum field solution 
{—0.5, 0.8} . Note the minimum of the energy surface at ci ~ 0.15. 



£ (fi, ci) in the parameter ranges [s, G {4, 9} Aym and ci G 
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FIG. 5: The energy density e (ci) = e (fi c (ci) , c\) at the phase transition point [i c in units of Aym- 



